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ABSTRACT 

I measure the physical properties of 38 transiting extrasolar planetary systems, bringing the 
total number studied within the Homogeneous Studies project to 82. Transit light curves are 
modelled using the JKTEBOP code, with careful attention paid to limb darkening, orbital ec- 
centricity and contaminating light. The physical properties of each system are obtained from 
the photometric parameters, published spectroscopic measurements and five sets of theoreti- 
cal stellar model predictions. Statistical errors are assessed using Monte Carlo and residual- 
Oh . permutation algorithms and propagated via a perturbation algorithm. Systematic errors are 

■ estimated from the interagreement between results calculated using five theoretical stellar 

(-H \ models. 

Qj The headline result is a major upward revision of the radius of the planet in the OGLE- 

- TR-56 system, from 1.23-1.38 Rjup to 1.734 ± 0.051 ± 0.029 Rjup (statistical and sys- 

tematic errors, respectively). Its density is three times lower than previously thought. This 
change comes from the first complete analysis of published high-quality photometry. Signifi- 
cantly larger planetary radii are also found for Kepler-15, KOI-428, WASP-13, WASP-14 and 
WASP-21 compared to previous work. 

I present the first results based on Kepler short-cadence data for Kepler-14, Kepler-15 
and KOI- 135. More extensive long-cadence data from the Kepler satellite is used to improve 
the measured properties of KOI-196, KOI-204, KOI-254, KOI-423 and KOI-428. The stellar 
component in the KOI-428 system is the largest known to host a transiting planet, at 2.48 ± 
0.17 ± 0.20 Rq. Detailed analyses are given for HAT-P-3, HAT-P-6, HAT-P-9, HAT-P-14 and 
WASP-12, based on more extensive datasets than considered in previous studies. 

Detailed analyses are also presented for the CoRoT systems 17, 18, 19, 20 and 23; Kepler- 
7, 12 and 17; KOI-254; OGLE-TR- 1 1 1, 1 13, 132 and L9; and TrES-4. 

I revisit the correlations between orbital period and surface gravity, and orbital period 
and mass of the transiting planets, finding both to be significant at the 4(T level. I conclude by 
discussing the opportunities for follow-up observations, the sky positions and the discovery 
rate of the known transiting planets. 

Key words: stars: planetary systems — stars: fundamental parameters 



1 INTRODUCTION 

Even though the first transiting extrasolar planet (TE P ) was 

discovered only as recently as 1999 jHenrv et alj llOOOl : 

ICharboimeau et all |2000|) . the number of recognised TEPs - 
and the size of the literature on them - has grown very quickly. 
The known planetary systems now display a bewildering diversity 
of qualities: from super-Earths with extremely short periods 
(CoRoT-7b; iLeger et alj 12009 ) to massi ve planets on wid e 
and highly eccentric orbits (HD 80606 b; iHebrard et all l2010h : 
a system of six planets showing strong mutu al gravitational 
perturbations (Kepler- 1 1 ; iLissauer et 'Zr i201lh ; circumbinary 
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planets toovle et al] l201ll ; IWelsh et al] l2012h ; and ones with a 
size inexplicably approaching twice that of J upiter (WASP- 17 b; 
lAnderson et al.ll20ld ; ISouthworth et alj2012bl) . 

Part of the reason for the known TEP population being such a 
menagerie is the variety of methods deployed in catching them. 
Many of the TEPs orbiting brighter stars (including the first, 
HD 209458 b) were found by radial velocity surveys and later 
shown t o transit. Oth ers (such as the second known TEP OGLE- 
TR-56 b; |Konacki eTa l. 2003) were unearthed via deep photometric 
searches and subsequently confirmed spectroscopically. The dom- 
inant population of TEPs currently stems from large-scale photo- 
metric surveys of stars of intermediate brightness (V = 9-13) cou- 
pled with spectroscopic follow- up programs, such as those of HAT 
dBakos et al.,,2002.) and WASP dl^ollacco et alj|2006h . This popula- 
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tion is strongly biased towards larger planets with orbital periods 
forb ^ lOd, as these Hot Jupiters show the greatest photometric 
variability. More recently, surveys such as CoRoT and Kepler have 
found success by obtaining uninterrupted photometry of stars us- 
ing space satellites. The objects found in this way - like the deep 
ground-based surveys - mostly orbit fainter stars, which hinders 
attempts to obtain follow-up observations to refine their physical 
properties or even spectroscopically confirm the planetary nature 
of the transiting body. 

An important problem with the study of TEPs is that their 
physical properties cannot be obtained simply by measuring some 
observable quantities and putting them through standard formulae: 
the number of quantities measurable directly from observations is 
one too few. An additional constraint is therefore needed, which is 
usually obtained by requiring the physical properties of the host star 
to match the predictions of stellar theory. In practise, there are sev- 
eral ways of implementing this constraint, and multiple sets of the- 
oretical predictions from which to choose. As different researchers 
select different approaches, their results are inhomogeneous and 
therefore non-trivial to compare. 

The diversity of the known TEPs, coupled with the varied 
methods used in their analysis, makes a consistent picture of their 
properties more difficult to attain. The current series of papers is an 
attempt to perform a uniform analysis for all suitable TEP systems, 
in order to produce homogeneous measurements of their physical 
properties. Such results are valuable in statistical investigations of 
the characteristics of this population (such as work similar to that 
bv lEnochet"ai]|2012l) . as well as providing an independent con- 
firmation (or otherwise) of published numbers. In some cases it 
is possible to achieve improved results simply by considering all 
available data rather than concentrating on only one set of observa- 
tions. 

In Paper I ( ISouthworthll2008h I presented the methods used 
for analysing the transit photometry of TEP systems, and applied 
them to the fo urteen object whic h had good light curves at that 
point. Paper II dSouthwortriboogl) discussed the more thorny is- 
sue of how to turn the measured photometric and spectroscopic 
parameters into physical properties; the primary problem here is 
the need to incorporate the predictions of theoretical stellar evolu- 
tionary models (or some other additional constraint) into the pro- 
cess. This was achieved for the fourteen TEPs by using three sets of 
stellar models, plus an empirical stellar mass-radius relation as an 
alternative constraint. In Paper III (Southworth 201G) the net was 
widened to include the full analysis of 30 TEPs with five sets of 
stellar models, plus a description of the methods used to incorporate 
conta minating light and orbital eccentricity into the analysis. Pa- 
per IV ('Southw ortlj201lh further extended the project to include 58 
TEPs, concentrating on those which had space-based light curves 
from CoRoT, Kepler, or the NASA Epoxi satellites. An improved 
empirical constraint was also used (based on stellar density, temper- 
ature and metal abundance), calibrated on the measured properties 
of 180 stars in 90 detached eclipsing binaries. 

In this work I present new photometric analyses of 30 more 
TEPs, and measure the physical properties of these plus a fur- 
ther eight objects previously considered in the current project. This 
brings the total number of TEPs within the Homogeneous Studies 
project to 82. Whilst this is by far the largest source of homoge- 
neous physical properties, it still represents less than half of the 
known TEP systems. The physical properties, follow-up status, sky 
positions and discovery rate of the known TEPs are finally dis- 
cussed in Sections [6] to [8] As with previous papers in this series. 



extensive tables of intermediate results are available for inspection 
in the online-only Supplementary Information. 



2 LIGHT CURVE ANALYSIS 

I have studied the light curves of each TEP system using the meth- 
ods discussed in detail in Paper I. The JKTEBOlf^code was used to 
find the best-fitting model for each light curve, from which the pho- 
tometric parameters were extracted. JKTEBOP represents the star 
and planet as biaxial spheroids with shapes governed by the mass 
ratio. The results in this work are all very insensitive to the adopted 
mass ratio, in part because the projected surfaces of the stars suffer 
the least distortions from sphericity close to the times of transit and 
occultation. The correct terminology for transits (primary eclipses) 
and occultations (secondary eclipses) was discussed in sect. 4 of 
Paper IV. 

The main parameters for JKTEBOP are the orbital inclination 
(i), and fractional radii of the star {ta) and planet (rb), defined as: 



TA = 



Ra 



r-b 



Rh 

a 



(1) 



where a is the orbital semimajor axis, and _Ra and _Rb are the 
volume-equivalent stellar and planetary radii. In this work I fitted 
for the sum and ratio of the fractional radii: 



rA + r-b 



rb 

TA 
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Ra 
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as these are less strongly correlated. In all cases the time of transit 
midpoint (To) was included as a fitted parameter. When modelling 
one transit, the orbital period (Porb) was fixed to a value taken from 
the literature. When modelling multiple transits at the same time, 
Porb was included as a fitted parameter. 

Limb darkening (LD) was considered in detail. Every light 
curve was solved with each of five different LD laws (see Paper I) 
and with three different approaches to the LD coefficients: (1) both 
fixed ('LD-fixed'); (2) the linear coefficients (ua) fitted and the 
nonlinear coefficients (va) fixed but perturbec0 by ±0.10 in the 
error analysis simulations ('LD-fit/fix'); (3) both coefficients fitted 
('LD-fitted'). In a few cases where the photometric data are of lim- 
ited quality the LD-fitted alternatives were not calculated. Initial 
or fixed values for the LD coefficients were obtained by bilinear 
interpolation to the correct TeS and \ogg within tabulated theo- 
retical prediction^. Once the best of the three options (LD-fixed, 
LD-fit/fix, LD-fitted) was determined, the final light curve parame- 
ters were taken to be the weighted mean of the relevant results for 
the four nonlinear laws. 

Errorbars for the par ameters were obtained f rom 1000 Monte 
Carlo (MC) simulations j Southworth et alj|2004l) . for each of the 
solutions with the adopted approach to LD. The largest of these 
was retained, and an additional contribution added to reflect any 
variation in the parameter values from the four nonlinear LD laws. 



^ JKTEBOP is written in FORTRAN?? and the source code is available at 
|http : / /www, astro ■ keele ■ ac ■ uk/jkt /codes /jktebop ■ html) 

The size of the perturbation was incorrectly stated to be ±0.05 in Pa- 
per IV, rather than the true valu e of ±0. 1 0. 1 thank the referee of a previous 
paper I Sout hworth et alj2 012a) for bringing this to my attention. 

Interpolation within multiple tables of theoretical LD coefficients was 
performed using the JKTLD code, written in FORTRAN?? and available at 
|http : / / www .astro ■ keele ■ ac ■ uk/ jkt / codes / jktld . html| 
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Finally, alternative errorbars were calculated using a residual per- 
mutation (RP) algorithm (as implemented in Paper I), which ac- 
counts for correlated noise. The RP errorbars were adopted if they 

were greater than the MC equivalents. 

Any contaminating 'third' light (e.g. lDaemgenetai]|2009l) . 
denote d by L3, was included a s a constraint following the proce- 
dure of'South worth et alj ilOVj) . A detailed investigation of the is- 
sue of third light was given in Paper III. Orbital eccentricity was 
treated similarly jSouthworth et alj|2009cl) . with constraints on the 
combinations e cos uj and e sin uj preferred over those for eccen- 
tricity (e) and periastron longitude (oj) directly (see Paper III). Nu- 
merical integration over long exposure times is important for sev- 
eral TEPs, most notably those observed in long cadence by the Ke- 
pler satellite, and was dealt with in JKTEBOP following the scheme 
given in Paper IV. 



3 CALCULATION OF PHYSICAL PROPERTIES 

The analysis of transit light curves gives the orbital ephemeris 
(Poih, To) and the photometric parameters (va, rh and i). For 
each system we also have the orbital velocity amplitude of the 
host star (Ka)- These measured parameters alone do not lead to 
a unique solution for the physical properties of the system, so an 
additional constraint has to be sought from elsewhere. The possible 
constraints are outlined belowQ 

In each case, the physical properties calculated were the mass, 
radius, surface gravity and density of the star (Ma, Ra, log gA, Pa) 
and of the planet (A/b, Rb, Qb, Pb), a surrogate for the planetary 
equilibrium temperature: 



1/2 



(3) 



(4) 



and the lSafronovl ( Il972h number: 

2 [VorbJ ^ \Rb) \Ma) ^ rb Ma 
The set of physical constants used in the process above were 
tabulated in sect. 3.3 of Paper I\|f|, and are ver y close to the val- 
ues pr oposed to the astronomical community bv lHarmanec & Prlal 
( l201lh . 



3.1 Additional constraint from theoretical stellar models 

The standard approach in obtaining the additional constraint is to 
use the predictions of theoretical stellar evolutionary models. These 
provide the mass and radius of the host star, guided by determina- 
tions of its effective temperature (Tefi), metal abundance ( ) 
and light-curve-derived density jSeager & Mallen-Ornelasll2003h . 
The full physical properties of the system are then straightforward 
to calculate. I adopted the procedure introduced in Paper II, which 



* A small fraction of TEP systems have direct distance measurements 
from parallaxes, or ste llar diameter measurements from interferometry (e.g. 



Ivon Braun et alj2012h . These observations can supply the crucial additional 
constraint, but are not possible for most TEP host stars due to their distance 
from the Earth. 

^ The radius of the Sun and of Jupiter were wrong by a factor of ten in 
sect. 3.3 of Paper IV, due to a typographical error The correct values are 
6.95508 X 10* Rq and 7.1492 x 10^ Rjup, respectively. This does not 
affect any of the calculations in Paper IV or other analyses by the current 
author 



begins with estimating the velocity amplitude of the planet (Kb)- 
The other measured quantities (Porb, Ka, ta, rb, i, e) were used 
to deter mine the ensuin g physical properties using standard formu- 
lae (e.g. lHilditchll20oil) . The value of Kb could then be iteratively 
refined to minimise the figure of merit 



fom 



(obs) 



^iA ) 



T 



(obs) 



(calc) 



/T-T(obs)\ 



(5) 



where ri"''^-' and T^^''^' are quantities determined from observa- 



' A 

;j(calc) 



tions. ii^"""' and T^g 



off 
(calc) 



were obtained for a given Kb value by 
interpolating in the tabulations of theoretical models for the cal- 
culated mass and observed [|^] of the star. The outcome of this 
procedure is a set of physical properties corresponding to the best 
agreement between observed and calculated quantities. The solu- 
tion control parameter, Kb, represents the entirety of the input 
from theoretical stellar models. Because the evolution of a star is 
highly non-linear with time, the above procedure was carried out 
for a range of ages starting from 0.01 Gyr and incrementing in 
0.01 Gyr chunks until the star was significantly evolved (surface 
gravity log g < 3.5). This led to a set of overall best physical prop- 
erties plus a model-dependent estimate of the age of the system. 

The statistical uncertainties on the input parameters were 
propagated using a perturbation analysis (Southworth et al. 200^ 
implemented in the JKTABSDIM code. This had the additional ben- 
efit of yielding a complete error budget for every output parame- 
ter, allowing identification of how best to improve our understand- 
ing of each TEP. The use of theoretical stellar models also incurs 
a systematic error due to our incomplete understanding of stellar 
physics. I therefore ran JKTABSDIM solutions for each of five sets 
of theoretical m odels (see Paper III) : aflref iClaretll200 4 l2005l 



2006, 2007), jPemargue et alj2004h . Temmo ( Pietrinfe mi et al 



2004), VRSS jVandenBerg et all |200^ and DSEP ( Pot ter etal 



2008). This yielded five sets of values and statistical errorbars for 
each output parameter. For the final result for each parameter I have 
taken the unweighted mean of the five values. The largest of the in- 
dividual statistical errors was adopted as the final statistical error, 
and the standard deviation of the five values was taken to represent 
the systematic error. 

Three of the derived physical properties have no depen- 
dence on stellar theory so are quoted without systematic error- 
bars: flb jSo uth worth et al. 2007); pA to a good approximation 
dSeager & Mallen-Ornelasll2003ll : and r^'q (Paper III). 



3.2 Additional constraint from eclipsing binary relations 

In Paper II I found a way to avoid the use of stellar models en- 
tirely, by defining an empirical mass-radius relation based on well- 
studied eclipsing binaries. This approach was not very successful 
because it did not allow for either the chemical composition or 
evolution of the host star. It was also subject to the known phe- 
nomenon that the components of low-mass eclipsing binaries tend 
to have lar ger radii and lower Tegs than those pred icted by stellar 
theory (e.g. lL6pez-Moralesll2007l : lRibas et al]|2008l) . 

An alternative approach is to explicitly include a dependence 
on Tcff and [^] . i Enoch et al.i (2010) proposed a set of calibra- 
tion equations giving log M or log 7? as functions of p. Tea and 
[^] , where p was chosen because a very good approximation 
to this quantity is calculable d irectly from transit light curves 
( Seager & Mallen- Ornelasl l2003'). In Paper IV I adopted this ap- 
proach, but with modified calibration coefficients obtained from a 
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Table 1. Final photometric parameters for the 38 TEPs analysed in this work. Orbital periods ai'e also given, and are either from this work or taken from the 
literature. Systems for which a constraint on orbital eccentricity or third Hght was required, are marked using a * in the column headed "e?" or "L3?". 



System 


Orbital period 

(days) 


e? 


L3? 


Orbital incHnation, 
i (degrees) 


Fractional stellar 
radius, T'a 


Fractional planetary 
radius, r\, 


Reference 


CoRoT-2 


1.7429935 (10) 


-*- 




87.45 ± 0.34 


0.1478 ±0.0023 


0.02462 ± 0.00035 


Paper IV 


CoRoT-17 


3.7681 (3) 






89.31°;^ 


n 1 c;7+0-043 


0100+° '"'^"' 

U.U1UU_Q QQQg 


This work 


CoRoT-18 


1.9000693 (28) 






86.8 ± 1.7 


0.1502 ± 0.0072 


0.0209 ±0.0013 


This work 


CoRoT-19 


3.89713 (2) 






sq 7+0-3 
' -3.6 


n 14-5+0-028 


01 n +0-0024 


This work 


CoRoT-20 


9.24285 (30) 






83.5 ±3.8 


0.070 ±0.019 


0.0062 ±0.0014 


This work 


CoRoT-23 


3.6313 (1) 






84 7+^-^ 

"^•'-3.0 


n 1fiS+0.022 




This work 


HAT-P-3 


2.8997360 (20) 






86.15 ±0.19 


0.1053 ±0.0020 


0.01 178 ±0.00030 


This work 


HAT-P-6 


3.8530030 (14) 






84.88 ±0.51 


0.1346 ± 0.0054 


0.01271 ±0.00070 


This work 


HAT-P-9 


3.922814(2) 






86. 10 ±0.54 


0.1 177 ±0.0060 


0.01250 ±0.00089 


This work 


HAT-P-14 


4.6276691 (50) 






83.01 ±0.27 


0.1211 ±0.0036 


0.00954 ± 0.00045 


This work 


Kepler-7 


4.8854922 (38) 




* 


85.32±0.16 


0.1491 ±0.0013 


0.01246 ±0.00015 


This work 


Kepler- 12 


4.4379637 (2) 




-k 


89.02+^;25 


n 19ASR + 0-00110 


n niAf;s7+o. 000071 

O.U14b87_|, QQQQgg 


This work 


Kepler- 14 


6.7901230 (43) 




-k 


86.64 ± 0.59 


0.1260 ± 0.0058 


0.00698 ±0.00029 


This work 


Kepler- 15 


4.942782(13) 




-k 


86. 144 ±0.051 


0.10002 ±0.00048 


0.010575 ±0.000086 


This work 


Kepler- 17 


1.4857108 (2) 




k 


89.73l»;2J 


n 1755+0-0018 


02403+° "°"2'' 

U.UZ1UO_Q QQQJ3 


This work 


KOI- 135 


3.0240933 (27) 




k 


84.89 ±0.21 


0.1409 ±0.0020 


0.01212 ± 0.00025 


This work 


KOI- 196 


1.8555588 (15) 




k 


89.5i?;^ 


0.1508]tQ^QQ20 




This work 


KOI-204 


3.2467220 (73) 




k 


84.30 ±0.83 


0.1485 ± 0.0090 


0.01259 ± 0.00098 


This work 


KOI-254 


2.45524122 (65) 


-*- 


k 


87.48 ±0.38 


0.0849 ± 0.0040 


0.01616 ± 0.00091 


This work 


KOI-423 


21.087168 (21) 


-*- 


k 


89.06 ±0.21 


0.0370 ±0.0015 


0.00339 ±0.00017 


This work 


KOI-428 


6.8731697 (96) 




k 


85.47 ±0.94 


0.1437 ± 0.0072 


0.00859 ± 0.00052 


This work 


OGLE-TR-56 


1.21191096 (65) 






73.47 ±0.37 


0.3292 ± 0.0048 


0.03378 ± 0.00059 


This work 


OGLE-TR-111 


4.0144477 (16) 






88.50 ±0.35 


0.0817 ±0.0018 


0.01033 ±0.00030 


This work 


OGLE-TR-113 


1.43247425 (34) 






87.80 ±0.86 


0.1572 ±0.0020 


0.02296 ± 0.00039 


This work 


OGLE-TR-132 


1.68986531 (64) 






83.83 ±0.96 


0.2056 ± 0.0078 


0.01941 ±0.00096 


This work 


OGLE-TR-L9 


2.48553417 (64) 






82.19 ±0.21 


0.1722 ±0.0025 


0.01928 ± 0.00028 


This work 


TrES-4 


3.5539268 (32) 




k 


82.33 ±0.47 


0.1697 ±0.0064 


0.01650 ± 0.00059 


This work 


WASP-1 


2.5199464 (8) 






88.0 ±2.0 


n 17Q7+0.0057 
u.±loi _0.0089 


0.0182lO;ri'I 


Paper 1 


WASP-2 


2.15222144(39) 




k 


84.81 ±0.17 


0.1238 ±0.0018 


0.01643 ±0.00030 


Southworth et al. (2010) 


WASP-4 


1 33823144 (32^ 






89.0 ± 1.0 


1825+'^°°" 


02812+° '""'22 


Southworth et al. (2009b) 


WASP-5 


1.6284246(13) 






85.8± 1.1 


0.1 847 ±0.0061 


0.02050 ±0.00091 


^uthworth et al. (2009a) 


WASP-7 


4.9546416 (35) 






87.03 ±0.93 


0.1 102 ±0.0061 


0.01053 ± 0.00070 


Southworth etal. (2011) 


WASP- 12 


1.0914222(11) 






83.3 ± 1.1 


0.3260 ± 0.0052 


0.03777 ± 0.00088 


This work 


WASP- 13 


4.353011 (13) 






84.88 ±0.45 


0.1382 ±0.0049 


0.01310 ± 0.00059 


This work 


WASP- 14 


2.2437661 (11) 






81.1±1.5 


0.2084 ± 0.0099 


0.0210 ±0.0010 


This work 


WASP- 18 


0.94145181 (44) 


* 




85.0 ±2.1 


0.2795 ± 0.0084 


0.0272 ±0.0012 


Southworth et al. (2009c) 


WASP-2 1 


4.3225060 (31) 






86.77 ±0.45 


0.1069 ±0.0037 


0.01 170 ±0.00054 


This work 


XO-2 


2.6158640 (16) 






88.8± 1.2 


n 19-37+0.0024 


oi^OO+O ''''''^^ 


Paper III 



larger and better-conditioned sample of eclipsi ng binary systemp 
A similar re-calibration has been performed bv^ Gillon et al.l ( l2()lll) . 

In the current work I have calculated the physical properties 
for each TEP system using equation 7 in Paper IV with the coef- 
ficients of the calibration for logjQ R covering stellar masses of 
0.2-3.0 Mq. These results are presented as an alternative to the de- 
fault solutions incorporating the predictions of theoretical stellar 
models. 

4 DATA ACQUISITION 

The CoRoT data used here are the N2 public version obtained in 
FITS format from the public archival]- The data for each object ex- 

^ See the catalogue of well-studied detached eclipsing binary star systems: 
|htt p : / /www . astro . ke ele ■ ac ■ uk/~ jkt/debcat/l 
[http ://idoc-corot.ias. u-psud. f r/| 



ist in short and/or long cadence flavours, with effective integration 
times of 32 s and 512 s respectively. I used the total flux as given by 
the "WHITEFLUX" column. The chromatic light curves were not 
considered as they add little to the results and have passbands which 
vary from star to star. I rejected any datapoints which were flagged 
as being unreliable (i.e. I required STATUS = 0). Numerical inte- 
gration was used for the long-cadence data, adopting A'^int = 3 as 
the number of integration points. 

The Kepler data used were the public data obtained from 
the Multimission Archive at STScI (MASlEl). These also come 
in two forms: short cadence (effective integration time 58. 84876 s; 
Gilliland et al. 2010) and long cadence (29.4244 min; Jenk ins et al.l 
20101) . Numerical integration is crucial for the interpretation of the 
latter type of data (see Paper IV for a demonstration) and I adopted 
A'int = 10 in all cases. The Kepler data cover Quarter (QO) to 

® |http : //archive ■ stsci ■ edu/kepler/data^earch/| 
search . php 
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Table 2. Spectroscopic quantities for the TEP host stars, taken from the hterature and adopted in the analysis in the present work. 



5 



System 



Velocity amplitude ( m s ^ ) 



Tcft (K) Reference 



Reference 



CoRoT-2 

CoRoT-17 

CoRoT-18 

CoRoT-19 

CoRoT-20 

CoRoT-23 

HAT-P-3 

HAT-P-6 

HAT-P-9 

HAT-P-14 

Kepler-7 

Kepler- 12 

Kepler- 14 

Kepler- 15 

Kepler- 17 

KOI- 135 

KOI- 196 

KOI-204 

KOI-254 

KOI-423 

KOI-428 

OGLE-TR-56 

OGLE-TR-111 

OGLE-TR-113 

OGLE-TR-132 
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,-1-4.3 
-4.4 
a+27 
-25 
7-1-8.5 
9.1 



603 ± 18 

312.4 ±29.0 
590 ± 14 
126 ±6 
454 ±9 
377 ± 34 

89.1 ±2.0 

115.5 ±4.2 
84.7 ±7.9 

218.9±5.8 
42.9 ±3.5 
48. 2^ 
683j 
78. 7 j 
399 ±9 
375 ±13 
85 ±11 

124 ±5 
110± 10 

125ll^« 
179 ±27 
212 ±22 
78± 14 
267 ± 34 
167 ± 18 
510± 170 
97.4 ±7.2 

125 ±5 

153.6 ±3.0 
242.1 ±3.0 

268.7 ± 1.8 
97± 13 

226 ±4 
55.7 ±5.5 
990 ±3 
1816.7 ± 1.9 
37.2± 1.1 

92.2 ± 1.7 



Gillonetal. (2010) 
Csizmadia et al. (201 1) 
Hebrardetak (2011b) 
Guenther et al. (2012) 
Deleuil et al. (2012) 
Rouan et al. (2012) 
Torres et al. (2007) 
Noves et al. (2008) 
Shporer et al. (2009a) 
Torres et al. (2010) 
Latham et al. (2010) 
Fortnevetal. (2011) 
Buchhave et al. (2011) 
Endlet al. (2011) 
Bonomo et al. (2012) 
Bonomo et al. (2012) 
Santerne et al. (201 la) 
Bonomo et al. (2012) 
Johnson et al. (2012) 
Bouch y et al, (2011) 
Sante rne et al. (2011b) 
Bouchy et al. (2005) 
Pont et al. (2004) 
Torres et al. (2008) 
Moutou et al. (2004) 
Snellen et al. (2009) 
Mandushev el^ (2007) 
Albrechtetak (2011) 
Triau d et al. (2010) 
Tria ud et al. (2010) 
Triau d et al. (2010) 
Helli er et al. (2009b) 
Hebb et al. (2009) 
Skillen et al. (2009) 
Blecic etak (2011) 
Triaud et al. (2010) 
Bo uchy etal. (2010) 
Naritaetal. (2011) 



5598 ±50 
5740 ± 80 
5440 ± 100 
6090 ±70 
5880 ±90 
5900 ±100 
5185±50 
6570 ±80 
6350 ±150 
6600 ± 90 
5933 ±50 
5947 ± 100 
6395 ±60 
5595 ± 120 
5781 ±85 
6041 ± 143 
5660 ± 100 
5757 ±134 
3820 ± 90 
6260 ± 140 
65 10 ±100 
61 19 ±62 
5044 ± 83 
4790 ± 75 
6210±59 
6933 ±58 
6200 ± 75 
6213±51 
5170 ±60 
5540 ±55 
5770 ±65 
6520 ±70 
6250 ± 100 
5826 ± 100 
6475 ± 100 
6455 ±70 
5800 ±100 
5340 ±50 



Schroteret al. (2011) 
Csizmadia et al. (201 1) 
Hebrardetak (2011b) 
Guenther et al. (2012) 
Deleuil et al. (2012) 
Rouan et al. (2012) 
Torres et al. (2007) 
Noves et al. (2008) 
Shporer et al. (2009a) 
Torres et al. (20W) 
Latham et al. ( 2010) 
Fortnevetal. (2011) 
Buchhave etak (2011) 
Endletak (2011) 
Bonomo et al. (2012) 
Bonomo et al. (2012) 
Santerne et al. (2011a) 
Bonomo et al. (2012) 
Johnson et al. (2012) 
Bouchv etak (2011) 
Santerne et al. (201 lb) 
Santos et al. ( 2006) 
Santos et al. (2006) 
Torres et al. (200^ 
Gillon et al. (2()07) 
Snellen e t al. (2009) 
Sozzetti et al. (2009) 
Albrecht et al. (2011) 
Sect. l53T] 

Max ted etak (2011) 
Maxtedetak (2011) 
Maxtedetak (2011) 
Fossati et al. (2010) 
Skillen et al. (2009) 
Joshi et al. (2009) 
Maxted etal. (2011) 
Bouchv et al. (2010) 
Burke et al. (2007) 



0.04 ± 0.05 

0.0 ±0.1 
-0.1 ±0.1 

-0.02 ±0.10 
0.14 ±0.12 
0.05 ±0.1 
0.27 ± 0.05 

-0.13 ±0.08 
0.12 ±0.20 
0.11 ±0.08 
0.11 ±0.05 
0.07 ± 0.04 
0.12 ±0.06 
0.36 ± 0.07 
0.26 ±0.10 
0.33 ±0.11 

-0.09 ±0.16 
0.26 ±0.10 
0.20 ±0.10 

-0.29 ±0.10 

0.25 ±0.08 
0.19 ±0.07 
0.09 ±0.08 
0.37 ±0.07 

-0.05 ± 0.20 
0.14 ±0.09 
0.17 ±0.05 
0.04 ± 0.05 

-0.03 ±0.09 
0.09 ±0.09 
0.00 ±0.10 
0.32 ±0.12 
0.0 ±0.2 
0.0 ±0.2 
0.00 ± 0.09 

-0.46 ±0.11 
0.45 ± 0.05 



Schroter et al. (201 k) 
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Hebrard et al. (201 lb) 
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Deleuil et al. (2012) 
Rouan et al. (2012) 
Torres et al. (2007) 
Noves et al. (2008) 
Shporer et al. (2009a) 
Torres et al. (2010) 
Latham et al. (2010) 
Fortn ev et al. (2011) 
Buchhave etak (2011) 
Endlet al. (2011) 
Bonomo et al. (2012) 
Bonomo et al. (2012) 
Santerne et al. (2011a) 
Bonomo et al. (2012) 
Johns on et al. (2012) 
Bouchv etak (2011) 
Santerne etal. (2011b) 
Santos et al. (2006^ 
Santos e t al. (2006) 
Torres et al. (2008) 
Gillon et al. (2007) 
Snellen et al. (2009) 
Sozzett i et al. (2009 ) 
Albre cht et al. (2011) 
Albrecht et al. (2011) 
Gillon et al. (2009) 
Gillon et al. (2009) 
HelUer et al. (2009b) 
Fossati et al. (2010) 
Skillen et al. (2009) 
Joshi et al. (2009) 
Hellier et al. (2009a) 
Bouchv et al. (2010) 
Burke et al. (2007) 



Quarter 2 (Q2) with additional data up to Quarter 6 (Q6) available 
for some objects. 

In order to account for slow variations in these data, both as- 
trophysical and instrumental, all data more than one transit duration 
away from a transit were discarded. A straight line was then fitted 
to the data around each transit (not including the datapoints in tran- 
sit) and divided out to normalise the transit to unit flux. A straight 
line was in almost all cases sufficient to remove the slow varia- 
tions, and is also unable to distort the shape of the transit. A few 
instances where a polynomial was used instead are noted below. In 
those cases where a large number of datapoints remained, they were 
converted into orbital phase, sorted, and then phase-binned into a 
much smaller number of normal points. The amount of binning was 
carefully chosen to avoid smearing out the transit shape. Details of 
this process are given below when relevant. 

Other data were obtained from public archives, from pub- 
lished papers, and by personal communication when necessary. 



5 RESULTS FOR INDIVIDUAL SYSTEMS 

In this section I present JKTEBOP and JKTABSDIM analyses of 38 
TEPs based on published high-quality data. The results were ob- 
tained using the same methods as those in previous papers, leading 
to homogeneous measurements for a sample of 82 TEPs. The final 
photometric parameters of the 38 TEPs considered in the current 
work are aggregated in Table[T] and their adopted spectroscopic pa- 
rameters are given in Table|2] I have enforced minimum errorbars 
of ±50 K for Tcft and ±0.05 dex for [^] , as the stellar effective 
temperature and metallicity scales set a lower limit on how pre- 
cisely we can measure these q uantities. Such minimum errorbars 
may still be optimistic (see e.g. lSruntt et alJ|2010aLl201^ . and this 
will be revisited in future. 

5.1 CoRoT-17 

Discovered bv lCsizmadia et aki i201ll) . the CoRoT-17 system har- 
bours a fairly massive planet (2.46 ± 0.25 Mjup) around an old 
star (S.Ola g Gyr). The only available light curve is that from the 
CoRoT satellite, which is comparatively poor due to the faintness of 



© 0000 RAS, MNRAS 000. 000-000 



6 John Southworth 



1 .00 



0.99 



o 0.98 ■ . .* ■' " • -. .. . • ■ . ^ 

' 0.97 ^^^^^^^^i 

0.96 F I , , , I , , , I , . , I , , , I . , , I , , , I - 



-0.06 -0.04 



-0.02 0.00 0.02 
Orbital phase 



0.04 



0.06 



Figure 1. Phased CoRoT 512 s cadence light curve of CoRoT-17 compared 
to the best fit found using JKTEBOP and the quadratic LD law. The resid- 
uals are plotted at the base of the figure, offset from unity. The purple line 
through the data shows the best-fitting model without numerical integration, 
and through the residuals shows the difference between this model with and 
without numerical integration. In this case the purple lines are barely visible 
behind the other lines. 
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Figure 2. The CoRoT and Euler telescope light curves of CoRoT-18 com- 
pared to the JKTEBOP best fits. Other comments are the same as Fig.[T] 



the star (V = 15.46). Icsizmadia et al] ( l201lh found that their ob- 
servations were consistent with a circular orbit, and that the CoRoT 
light curve is contaminated with a third light of L3 = 8 ± 4 %. 

The CoRoT data comprise 10 248 512 s cadence datapoints 
covering 23 transit events. 1571 datapoints in the vicinity of transits 
were retained (see Sect.|4l(, of which ten were subsequently rejected 
by a 4(j chp. I modelled the light curve using numerical integration 
(Mnt = 3) and with Po^b and To included as fitted parameters. 
LD-fixed and LD-fit/fix solutions were attempted, but LD-fitted so- 
lutions were avoided due to the modest quality of the light curve. 
The Monte Carlo algorithm implemented in JKTEBOP includes the 
application of random perturbations to the starting parameters for 
each simulation, in order to sample the parameter space well. I had 
to turn this option off as the perturbations yielded many poorly con- 
ditioned results, a situation also encountered for Kepler-4 in Pa- 
per IV (sect. 6.2). The RP errorbars were therefore adopted, as they 
are significantly larger than the MC errorbars. The full light curve 
solutions are given in Table Al. 

A summary of the light curve solutions is given in Table A2 
and the best fit to t he CoRoT light curve is plotted in Fig.[T] The 
parameters found bv lCsizmadia et af] ( l201lh agree with my own to 
within the errorbars, although I find a solution with slightly larger i, 
smaller ta and rb, and asymmetric errorbars. The physical proper- 
ties of CoRoT- 17 are given iri Table A3 an d again show good agree- 
ment with those of Csiz madia^et alj ( 1201 ih . except that the stellar 
surface gravity found in that work (log g = 4.40 ± 0.10) is incom- 
patible with the star's mass and radius (Ma = 1.04 ± 0.10 
and Ra = 1.59 ± 0.07 Rq). I provide the first measurement of the 
surface gravity of the planet. 

CoRoT-17 would benefit from new photometric and spectro- 
scopic observations, in order to refine its physical properties and 
check the spectroscopic T^s value. The orbital ephemeris is impre- 
cise, and will be unable to accurately predict transit midpoints by 
2012 August (see Sect.|7ll. 



5.2 CoRoT-18 

The d iscovery of CoRoT-18b was presented by iHebrard et al.l 
Ih) . and contains a relatively massive planet (Mb = 3.27 ± 
0.17 Mjup) orbiting an active star with a rotation period of 5.4 ± 
0.4 d. Its relative faintness (V = 15.0) made RV measurements 
difficult. This was not a problem for the Ka value, which is large , 
but observations o f the Rossiter-McLaughlin effect teossitetll9M 
|m cLaughlinl 1 924) only constrain the orbit to be prograde and prob- 
ably aligned. 

The CoRoT data are all of 32 s cadence, and cover 13 consec- 
utive transits. Of the 56 819 original datapoints, 6932 in the vicinity 
of a transit were retained a nd then phase-binned b y a factor of ten 
to yield 694 normal points. iHebrard et alj ( l2011bl ) modelled these 
data after taking the unusual approach of binning them in time by 
a factor of 16 to be equivalent to 512 s cadence, to improve compu- 
tational tractibility. They therefore substantially lower the informa- 
tion content of the data they analysed. By comparison, the phase- 
binned data I worked on have a cadence equivalent to 32.2 s. These 
were modelled with the assum ptions of L3 = 2.0 ± 0.1 % and 
e = (Hebrar detal]|2011bh . The MC errors were found to be 
larg er than the RP ones, an d the full solutions are in Table A4. 

IHebrard etalj ( l201 id) also presented ground-based transit ob- 
servations of CoRoT- 1 8, from the Swiss Euler Telescope and Euler- 
Cam. This shows intriguing evidence around phase —0.008 for the 
passage of the planet over a starspot, but as noted by Hebrardetal] 
( l2011bh the deviation cannot conclusively be assigned to this phe- 
nomenon. It was conservatively treated as correlated noise in my 
analysis. G. Hebrard (private communication) confirmed that these 
data were obtained through an R filter. They were modelled as with 
the CoRoT light curve, except that third light was not included (Ta- 
ble A5). 

The light curve solutions are summarised in Table A6 and 
show that the CoRoT and Euler data agree well. The values from 
iHebrard et al . (201 lb) differ by about la, probably due to their ap- 
proach to binning the data. The best fits are shown in Fig.|2] 

The physical properties of CoRoT-18 are sho wn in Table A7 
and sho w that I find a smaller system scale than IHebrard et al.l 
feOUb l, meaning that the masses and r adii of both components 
are slightly smaller. iHebrard et aljj201 Ibh were unable to constrain 
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Figure 3. The CoRoT light curve of CoRoT-19 iGuenther et al.l2012l) com- 
pared to the JKTEBOP best fit. Other comments are the same as Fig.[T] 
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Figure 4. The CoRoT fight curve of CoRoT-20 compared to the JKTEBOP 
best fit. Other comments are the same as Fig.[T] 



the age of the star, as different age indicators were inconsistent with 
each other. The more precise hght curve parameters derived in the 
current work point to a large age of roughly 10 Gyr, which disagrees 
with the comparatively fast rotation and strong lithium abundance 
indicative of an age of < 1 Gyr. The evolutionary status of the star 
is therefore uncertain. 1 provide the first measurement of gh- The 
system would benefit from a better light curve, which would allow 
refined physical properties to be obtained and also shed light on the 
star's age via a more precise measurement of its density. 

5.3 CoRoT-19 

The planetary natur e of this system was discovered by 
iGuenther etail ( l2012h . The CoRoT data comprise 56 837 data- 
points at 32 s cadence covering seven consecutive transits, of which 
9229 points near transit were binned by a fact or of ten to give 924 
normal points (Fig.[3ll. iGuenther et alj ( 1201 2[) constrained L3 < 
0.32 % and e < 0.15 (3a) so third light and orbital eccentricity 
were neglected in my analysis (Table A8). The RP errors were used, 
as they are modestly larger than the MC ones. From Table A9 it is 
clear tha t I find a slightly diffe rent solution to the light curve com- 
pared to IGuenther et al.l ( I2OI2I) , with much larger and asymmetric 
errorbars. The two solutions are compati ble within my errors, but 
not within those of'Guent her et al.l ( I2OI2I) . 

The physical properties I find for the CoRoT- 19 system (Ta- 
ble A 1 0) are imprecise and agree with those of IGuenther et al.l 
( I2OI2I) . A discrepant solution with a young stellar age (zero-age 
main sequence versus ~5 Gyr) was found when using the mod- 
els, but restricting the possible solutions to older ages yielded re- 
sults which were consistent with those from the other four model 
sets. I provide the first measurement of g\^. Further photometry and 
spectroscopy of CoRoT-19 are needed to refine measurements of 
its physical properties, including the orbital ephemeris. 

5.4 CoRoT-20 

As found bv lDeleuiletd] ( l2012l ). CoRoT-20 contains a rather mas- 
sive planet (5.06 ±0.34 Mjup) in an eccentric (e = 0.562 ±0.013) 
and long-period (9.24 d) orbit around a G2V star. Third light is 
constrained to be less than 0.6% so was neglected in my analysis. 
I adopted e cos oj = 0.312 ± 0.022 and e sinoj = 0.468 ± 0.017 
jPeleuil et all2012l) . The CoRoT data comprise 56 855 datapoints 
at 32 s cadence but cover only three transits. The 2479 datapoints 



near these transits were binned by a factor of five into 496 normal 
points and modelled using JKTEBOP (Table All). The light curve 
is of li mited quality (Fig.B so I obtained imprecise results. Com- 
pared to lPeleuil et alj ( 120 12i) I find a solution with Icr lower i and 
correspondingly larger rA and r^, with much larger errorbars (Ta- 
ble A12). 

A similar problem was found, when calculating the physical 
properties of the system, as for CoRoT-19 above. The mod- 
els gave a discrepant solution with a very young age (Table A13). 
This was solved in the same way, by ar tifically restricting th e pos- 
sible ages for the Y^ models to > 1 Gvr IPeleuil et alj i2012h opted 
for the low-age solution, which is more consistent with their dif- 
ferent photometric parameters and with their detection of the Li 
6708 A absorption line. The discrepancy supports the much larger 
errorbars I find, and suggests further observations of CoRoT-20 
would be valuable. I find significantly different ph ysical proper- 
ties fo r the planet: Mb = 5.06 ± 0.36 ± 0.04 Mjup jPeleuil et all 
I2OI2I obtained 4.24 ± 0.23 M.T„p). i^h = 1.16 ± 0.26 Rjup 
(0.84 ± 0.04 Rjup) and T^'q = 1100 ± 150 K (1002 ± 24 K). I 
provide the first published measurement of g^- 

5.5 CoRoT-23 

CoRoT-23 dRouan et alj2012l) contains another massive and eccen- 
tric planet (Mb = 3.06±0.31 Mjup ande = 0.16±0.02), orbiting 
a somewhat evolved star (loggA ~ 4.011q5i)- CoRoT observed 
it at 512 s cadence initially (3598 datapoints covering nine transits) 
and then switched to 32s cadence (130544 datapoints covering 15 
transits) after it was identified as a planet candidate. The 584 dat- 
apoints at 512 s cadence in the region of transits were modelled 
using TVint ~ 3. The 20614 datapoints at 32 s cadence near tran- 
sit were phased and binned by a factor of 40 to yield 516 normal 
points. Both sets of data were modelle d adopting e — 0.16 ± 0.02, 
a; = 52° ± 9°, and with L3 = 7.2%. lRouan et al.l ( |2012|) did not 
give an errorbar for L3, so I used an estimate of 3%. For the 512 s 
data (Table A14) the RP errorbar was larger than the MC one for 
i, whereas for the 32 s data (Table A15) the RP errorbars were all 
smaller than the MC equivalents. 

The light curve solutions are summar ised in Table A16 a nd 
the best fits plotted in Fig.jS] Compared to iRouan et al.l ( l2012h . I 
find larger errorbars, and 1.5a larger rA and rb. The solutions 
to both datasets are ill-defined and have asymmetric errorbars. 
Pespite these issues, my physical properties of CoRoT-23 are in 
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Figure 5. The CoRoT light curve of CoRoT-23 compared to the JKTEBOP 
best fit. Other comments are the same as Fig.[T] 



good agreement with those of iRouan et alj j2012h both for val- 
ues and errorbars (Table A17). A lower Tctt would match the pho- 
tometric parameters better, and an improved light curve is im- 
portant for refining these parameters plus improving the orbital 
ephemeris. CoRoT-23 would be an interesting system to study fur- 
ther, due to its evolved host star, but its faintness may be an issue 
(V — 15.63 ± 0.07). I provide the first published measurement of 
5b. 
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Figure 6. Light curves of HAT-P-3 compared to the best fit found using 
JKTEBOP and the quadratic LD law. The residuals are plotted at the base 
of the figure, offset from unity. The light curve sources and passbands are 
labelled. 



5.6 HAT-P-3 

Presented bv lXorres et al] ( l2007h . HAT-P-3 is a fairly typical TEP 
system except that it is metal-rich ([§] = 0.27 ± 0.05) and 
the plane t is unexpe ctedly small, two qualities which may be con- 
nected. JbreeTeTn] I2OO7) found that a heavy-element core of mass 
~75 Mffi was needed for theoretical models to match the planet's 
radius-. IChan et alj ( 1201 lb s uggested ~ 100 Mq . 

The discovery paper jTorres et al.ll2007h included follow-up 
photometry of one transit obtained with the KeplerCam imager on 
the FLWO 1.2 m telescope. The dataset comprises 386 datapoints 
obtained in the i-band. I multiplied the errorbars of the data by 
Vs. 503 to obtain a reduced of Xv = 1-0. The data were mod- 
elled with JKTEBOP, with the findings that the LD-fit/fix results are 
best, RP errors are similar to MC errors, and the errorbars are sig- 
nific antly asymmetric (T able A 18). 

iGibson et alj ( |2008|) published an extensive analysis of HAT- 
P-3 based on photometry of seven transits obtained with LT/RISE, 
with the intention of probing for transit- timing variations (TTVs). 
To solve these data 1 phase-binned them by a factor of 40 to reduce 
the 19 681 original points to 493 normal points, with the rejection 
of 24 points by a 4(j clip, and multiplied the errorbars by \/2.358 
to obtain xl = 1-0. The LD-fit/fix solutions are best, MC and RP 
erro rs are similar, and the erro rbars are symmetric (Table A 19). 

iNascimbeni et al. I( l201lh obtained high-speed photometry of 
one transit of HAT-P-3 (exposure time 2 s, cadence 5 s) using 
AFOSC on the Asiago 1.8 m telescope. The 2247 original data- 
points were already condensed to 122 binned points in the datafile 
lodged with the CDS; I modelled these after multiplying the errors 



by Vl.985. The results (Table A20) show insignificant correlated 
noise and th at the LD-fit/fix so lutions are best. 

Finallv. lChariet observed HAT-P-3, again using Ke- 

plerCam, over three transits in the i-band and three more in the 2- 
band. The two sets of data (1258 and 506 points respectively) were 
solved separately with fixed Porb . I found that the RP errors were 
similar to the MC ones and that the LD-fit/fix solutions are best 
(Tables A21 and A22). The scatter of the z-band data is around 
1.3mmag, usefully lower than the 2.2mmag for the i-band data. 
The best fits to all datasets are shown in Fig.|6] 

The light curve solutions are summarised in Table A23, and 
show a comparatively poor agreeme nt between dataset s. Ths is pri- 
marily due t o the z-ba nd data from ' Chan et al] j201 ll) but also to 
the data from lTorres e t al. (2007). The Xv of the model "the values 
from different light curves agree" is 1.7-2.1 for rA -I- ^b, fc, ta and 
rb, and 1.0 for i. If the z-band data are discounted then I instead 
find a much better x^ = 0.3-1.0. The source of this discrepancy 
is not clear, so I have adopted the results ignoring the z-band data. 
The agreement between my final photometric parameters and those 
of previous studies is not great; my parameter values mostly lie in 
the middle of the distribution of those from previous studies. This 
situation is understandable in that I have included multiple datasets 
which are in modest disagreement, whereas previous researchers 
have tended to concentrate on only their own data rather than the 
full set of available observations. 

The results from my JKTABSDIM analysis (Table A24) show 
good agreement with previous work, except that I find modestly 
larger radii for both star and planet. The larger stellar radius shifts 
the age of the system from 0.4— 1.6Gyr to 7.5^3 ^ Gyr. My new 
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Figure 7. Light curves of HAT-P-6 compared to the JKTEBOP best fits. 
From top to bottom the datasets are 2:-band dNoves et al. 20Q8i). J-ban d 
)Todorovet"ani2012l) and J-band jTodorov et alJl2012l ; ISada et alj|2012h . 
Other comments are the same as Fig.|6] 



solution should be the most reliable as I have based it on all avail- 
able photometry (albeit with one dataset rejected). Further photom- 
etry would be useful in helping to understand the two discrepant 
datasets, but HAT-P-3 is already well-characterised so such effort 
is better spent on other TEPs. 

5.7 HAT-P-6 



HAT-P-6 b was discovered by iNoves et alj ( |2008[) and is known 
to have a retro grade orbit from Rossiter-McLaughlin observations 
jHebrard et 3111201 la), in agreement with the high T^ft of the host 



Todorov et alj ( I2OI2I) have presented Spitzer observations of 



star. 

occultations at 3.6 and 4.5 fj,m which provide a useful constraint on 
the orbital shape: e cos oi < 0.004 (3cr). 

Two z-band fo llow-up transits of HAT-P-6 were presented by 
iNoves et al] ( l2008h . observed with the FLWO 1.2 m and Kepler- 
Cam on two nights separated by 312 d. I used the decorrelated ver- 
sion of the data, as is plotted in fig. lb of that paper The number of 
datapoints is 499 and 515, respectively, and 1 multiplied the error- 
bars by V3.76 to obtain Xv = 1-0. The results (Table A25) show 
that the LD-fit/fix solutions are good and that the MC errors are 
larger than the RP alternatives. 

Aside from their Spitzer photometry, iTodorov et al ] ( l2012h 
presented two transits, one observed in the 7-band using a 0.36 m 
telescope and one observed in the J-band using the KPNO 2. 1 m 
telescope with FLAMINGOES infrare d imager and spec trograph. 
The latter dataset was also presented bv lSada et al.l ( I2OI2I) . who re- 
quired a fifth-order polynomial to normalise the light curve to unit 
flux due to the increased complications involved in infrared versus 
optical observations. 

The i-band light curve comprises 649 datapoints which suffer 
from a large scatter, requiring the use of fixed LD coefficients as 
well as a fixed out-of-transit brightness in order to arrive at a 'rea- 
sonable' model for the data (Table A26). For this reason I did not 
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Figure 8. The light curves of HAT-P-9 compared to the JKTEBOP best fits. 
The top three datasets are from Shp orer et all j2009bl) and the fourth is from 
iDittmann et alj j2012h . Other comments are the same as Fig.|6] 



use the results from these data when calculating the final photomet- 
ric parameters. 

For the 560 J-band datapoints, I multiplied the errorbars by 
\/2.65 and calculated JKTEBOP solutions using a fifth-order poly- 
nomial flux normalisation (Table A27). The RP errors are twice as 
big as the MC errors, indicating that these data contain significant 
correlated noise (see in particular around phase —0.04 in Fig.|7] I 
adopted the LD-fit/fix results and RP errors. The best-fitting LD is 
weaker than theoretically predicted, but only by la. 

The 2; -band and J-band results show a reasonable agreement 
(Table A28), except for k where = 2.1. This is likely due to the 
lower reliability of the J-band transit, and the errorbars of the fi- 
nal photometric parameters have been infl ated accordingly. Fu rther 
photometry of HA T-P-6 would be useful. iNoves et alj ( I2008h and 
[Torres et al. 1(12008') find almost identical photometric parameters to 
each other, which are also in agreement with my own. 

I find physical properties for the HAT-P-6 system which are 
in good ageement with previous work (Table A29) but indicate a 
sightly larger star and smaller planet. My new results are based on 
more extensive data than other studies of HAT-P-6 so should be the 
most reliable. Further improvement could be achieved by obtaining 
new photometry, as the light curves are the current bottleneck in 
our understanding of this TEP system. 



5.8 HAT-P-9 



The di scovery of HAT-P-9 b was announced by ISh porer et al.| 
( l2009bl) . who also obtained two follow-up transits with the FLWO 
1.2 m (2;-band), one with the WISE 0.4 6 m (unfiltered) and on e with 
the WISE 1.0 m (.R-band) telescopes. iMoutou et ahl ( 1201 ll) mea- 
sured the Rossiter-McLaughlin effect to find A = — 16±8, which is 
consistent with alignment between the orbital axis of the planet and 



© 0000 RAS, MNRAS 000. 000-000 



10 John Southworth 



the spin axis of tlie star. IPittmann et alj j2012l) obtained two more 
transit observations in the /-band using the 1.6 m Kuiper telescope. 

For the z-band data, 675 points over two transits, I fitted for 
Porb- The RP errors were smaller than the MC errors and full so- 
lutions are in Table A30. For the unfiltered (399 datapoints) and 
_R-band (154) data I did not fit for Porb- 

For the unfiltered data 1 adopted LD coefficients appropriate 
for the CoRoT white light curves as these are close enough for the 
current purposes, and found that the RP errors were slightly bigger 
than the MC ones (Table A3 1 ). For the P-band data I multiplied the 
errorbars by \/5.45 to get xi ~ 1-0, and found that the data prefer 
a central transit (Table A32) in disagreement to the other datasets 
for this star. The different light curve shape can clearly be seen in 
Fig.i 

The two /-band transits appear to have slightly different 
depths (see fig. 2 in Dittmann et al. 2012) but individual solutions 
showed that the difference is not significant. I therefore modelled 
all 1 144 datapoints together, with Porb included as a fitted quantity 
(Table A33). 

The light curve solutions are summarised in Table A34. 
Agreement between the three datasets (z, R, I) is poor (Xv ~ 2.7 
for ta + rh, Xv = 8.7 for fc) and not easily explicable. The R- 
band results appear to be somewhat different to the z and / ones, 
but rejection of these does not improve matters (xt ~ 13.7 for fc). 
I therefore calculated the final photometric parameters from the z, 
R and / individual results, with the errors inflated to account for 
the disagreement. Further photometry of HAT-P-9 is recommended 
to help understand this discrepan cy and improve the me asured sys- 
tem parameters. The results from lShporer et al.l ( l2009bl) agree with 
my values to within my errorbars, but not to within their errorbars. 

The JKTABSDIM results (Table A35) show a very good agree- 
ment between the values obtained via di f ferent s tellar models, and 
with the values found by IShporer et al] ( l2009b l). HAT-P-9 would 
benefit from further photometric and spectroscopic (Toff and [§] ) 
study. 

5.9 HAT-P-14 

Discovered by [Torres et al. U2O10h . the HAT-P-14 system contains 
a relatively massive planet (2.27 ± 0.08 Mjup). Its low orbital in- 
clination (i.e. high impact parameter) is potentially helpful in mea- 
suring the light curve parameters to high precision, and in probing 
for variations in thes e quantities. It was in dependently discovered 
by the WASP survey jSimpson et al.ll20lTi) . from whom additional 
photometry and spectroscopy is available. The spectroscopic prop- 
erties from the two consortia {Tea and [^^] ) agree. One more tran- 
sit lig ht cur ve has been presented by Nascimbeni et al. (201 1), and 
Iwinn et al.l ( l201lh have found it to be in a retrograde orbit from 
observing its Rossiter-McLau ghli n effect. 

Both lTorres et al.l ilOld) and lSimpson et alj ( 1201 ll) find a sig- 
nificant orbital eccentricity. This is expected given the correlation 
between t he presence of a relatively massive planet and an eccen- 
tric orbit dSouthwortheTa l. 2009c). I adopted e = 0.107 ± 0.013 
andw = 94° ± 4° dTorres et aU201ol : IWinn et al.ll201ll : |Pont erd] 
l201lh in the analysis below. Due to the limited quality of individual 
light curves I did not att empt LD-fitted solutions. 

[Torres et"al] j2010l) presented i-band observations of five tran- 
sits using the FLWO 1.2 m, for which I multiplied the errorbars by 
V9.82 to get xS = 1-0. The full solutions are in Table A36 and the 
MC errors are larg er than the RP ones. 

The work by Simpson et includes two V+R band 

transits obtained using LT/RISE, of which one is only patially cov- 
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Figure 9. Light curves of HAT-P-14 compar ed to the JKT EBOP best fit. 
From top to bottom the datasets are t-band fr om lTorres e t al 1 201(?) , FTN Z 
and R ISE V+R from I Sim pson et al.ll20l"lh and R from .Nascimbeni et alj 
i201lh . Other comments are the same as Fig.[6] 

ered, and one Z-band transit using FTN/Spectral. I put the for- 
mer onto a common flux zeropoint and multiplied the errorbars by 
V5.15 before the JKTEBOP analysis. For both datasets (Table A37 
and Table A38), the LD-fit/fix solutions returned unphysical LD co- 
efficients so I adopted the LD-fixed results. The RP errors are sig- 
nificantly larger than the MC errors, showing that correlated noise 

is important for both datasets. 

The data of lNascimbeni et si\ ( l201ll) were originally 2247 dat- 
apoints but the file deposited at the CDS is binned to 109 nor- 
mal points. Solutions of these data (Table A39) once more yielded 
unphysical LD coefficients so the LD-fixed solutions had to be 
adopted. 

A summary of the light curve solutions can be found in Ta- 
ble A40 and a plot of the best fits in Fig.|9] Agreement between the 
four light curves is within the uncertainties for va + '"b, i and ta, 
and slightly greater than the uncertainties for k and Vh- Literature 
value agree with my own, but the errorbars are in some cases too 
small. This situation persists for the physical properties of the sys- 
tem (Table A41). More photometry and RV measurements would 
be useful, especially in reducing the reliance on theoretical LD co- 
efficients, but HAT-P-14 is a well-understood system. 



5.10 Kepler-7 

Kepler-7 contains a low-density TEP (0.094 ± 0.009 pjup) and a 
slightly evolved sta r (logg = 3.97 ± 0.02). The discovery paper 
dLatham et al.l2010l) presented an analysis of the Kepler QO and Ql 
data. The studies o f lKipping & BakosI ( 1201 ll) and Paper IV had ac- 
cess also to the Q2 data, although n one of the three datasets were 
short-cadence. iDemorv et al.l ( 1201 1[) presented a study of the QO- 
Q4 data, of which Q3 and Q4 are short-cadence, and found the 



© 0000 RAS, MNRAS 000, 000-000 



Studies of transiting extrasolar planets. V. 11 



1.000 




0.980 



-0,04 



-0.02 0.00 

Orbital phose 



0.02 



0.04 



Figure 10. The Kepler light curves of Kepler-7 compai'ed to the JKTEBOP 
best fits. The top light curve is the Q0-Q3 long-cadence data and the bottom 
one is the Q3 short-cadence data. Other comments are the same as Fig.[T] 



planet to have a high albedo from consideration of the occultation 
visible in the Kepler data. Their analysis included asteroseismolog- 
ical constraints, which are not used here in order to avoid inhomo- 
geneity. 

In the current work I have analysed the Q3 data, which are 
the only short-cadence data publicly available for Kepler-7. The 
improved time sampling of these data means that my results are 
better than those from Paper IV, which were based on long-cadence 
data. A third ligh t value of L3 = 1.4 % is listed in the Kepler Input 
Catalogue fKIC: Isrown et al1l20I ih . but without uncertainties. I 
therefore adopted an uncertainty of 0.5%. 

A new orbital ephemeris was derived from the Q0-Q3 long- 
cadence data, fitting only those data near one of the 42 transits and 
using Nint ~ 10. This yielded: 

To = BJD(TDB) 2 455 074.756840(53) + 4.8854922(38) x E 

where the reference time of minimum has been carefully chosen 
to be near the midpoint of the data interval, thus minimising its 
covariance with the orbital period. E represents an integer cycle 
count, and the bracketed quantities give the uncertainty in the last 
digit of the preceding number. Errorbars were calculated using MC 
and RP simulations and the larger of the two options is quoted. The 
two sets of errorbars were very similar, and also close to the formal 
errors calculated by JKTEBOP from the covariance matrix of the 
best-fitting solution. 

The Q3 data comprise 125 082 datapoints covering 17 full 
transits. One transit had a brightening feature affecting six dat- 
apoints around time BJD(TDB) 55133.35 which could indicate 
a stellar flare or instrumental artefact, so the data for this tran- 
sit were ignored. 16 897 datapoints around the remaining transits 
were phase-binned by a factor of 25 to give 676 normal points. 
These data were modelled using JKTEBOP (Table A42). The LD- 
fitted solutions are good enough to be adopted, and the RP errors 
are smaller than the MC ones (as usual for phase-binned data). 

The light curve solutions are summarised in Table A43 and 
the best fit to the data is plotted in Fig.[lO] The new results from the 
Q3 short-cadence data are in excellent agreement with those from 
Paper IV but are more precise. I combined the two to get the fi- 
nal photometric parameters. The values found bv lKipping & BakosI 
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Figure 11. The Kepler Ught curves of Kepler-12 compared to the JKTEBOP 
best fits. From top to bottom are the phase-binned short-cadence light curves 
from Q2 to Q6. Other comments are the same as Fig.|6] 



1 I2OI h are in agreement with my new results, whereas the ones 



Latham et alj ( l20IOh are in marginal disagreement 



from 

When determining the physical properties of Kepler-7 I found 
that the results using the different stellar models are in compara- 
tively poor agreement (Table A44), as often happens when the host 
star is evolved. My final results give a more massive star and planet 
than previously found, but within the errorbars. Kepler continues 
to observe this system, so even better results will available in the 
future. 



5.11 Kepler-12 

Discovered by iFortnev et al.l ( 1201 ih . the Kepler- 12 system con- 
tains the second-least-dense planet after WASP- 17 (* Anders on et al.l 
2010, 2011; Southworth et al. 2012b) and ahead of Kepler-7 (see 
above). iFortnev et alj ( I20I1 ) obtained the constraint e < 0.09 at 
the 3cr level. Data from QO to Q6 are publicly available, but the 
QO and Ql data are only long-cadence so I did not consider them 
further. Q4 is split into three periods, of which the Q4. 1 data are 
available but the Q4.2 and Q4.3 data are not. 

I analysed the Q2 to Q6 data, all of which are short-cadence, 
using the season-dependent third light values given in the KIC and 
assuming an uncertainty of 0.5%. The data from different quarters 
had to be analysed separately in order to be able to assign the ap- 
propriate third light values. In each case the transits were extracted 
from the full data, normalised to unit flux, and then phase-binned 
by a factor of 25 (a factor of seven for Q4 as it contains three times 
fewer datapoints). The solutions are given in Tables A45 to A49. In 
each case the MC and RP errors are similar. The light curve solu- 
tions are summarised in Table A50 and the best fits are plotted in 
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Figure 12. The Kepler Q3-Q6 short-cadence hght curve of Kepler- 14 com- 
pared to the JKTEBOP best fit. Other comments are the same as Fig.|6] 



Fig-El The five sets of solutions agree well with e ach other, and 
the co mbined results also agree well with those of iFortnev et alj 
(l201lh . 

The resulting physical pr operties (Table A51) are also in very 
good agreement with those of iFortnev et alj i201lh . Whilst Kepler 
continues to observe this system, the precision of its physical prop- 
erties is limited by the spectroscopic properties of the host star and 
by imperfect interpolation within tabulated predictions from theo- 
retical stellar models. 
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Figure 13. The Kepler Q3 to Q6 short-cadence light curves of Kepler- 15 
compared to the JKTEBOP best fits. Other comments are the same as Fig.|6] 



5.12 Kepler-14 

Discovered by iBuchhave et"aD bOllh . the system holds a mas- 
sive (7.7 ± 0.4 Mjup) hot Jupiter orbiting a comparatively hot star 
(Toff = 6395 ± 60 K). The Kepler light curve is strongly con- 
taminat ed by the light from a similarly bright star within 0.29" of 
system. ISuchhave et al.l ( l2oTlh measured a magnitude difference in 
the Kepler passband of Amxp = 0.44 ± 0.10, and corrected both 
the light curves and RVs for this contamination. Given the very low 
contamination levels from KIC (0.007-0.010 for the four Kepler 
seasons) it seems that the only relevant third light comes from the 
very close star. Its closeness means that the data from different sea- 
sons will be affected in the same way, so it is reasonable to adopt 
1/3 = 0.333 ± 0.062 when fitting all the Kepler data. 

Kepler observed Ke pler-14 in long c a dence in Q0-Q6, and in 
short cadence in Q3-Q6. iBuchhave et al.l ( 1201 ih opted to use the 
full set of long-cadence data only, for homogeneity and simplic- 
ity and because the transit is relatively long (L. Buchhave, private 
communication). Results based on short-cadence data stand a good 
chance providing improved results due to their much better time 
resolution. I have therefore modelled the Q3 to Q6 short-cadence 
data. 

These data comprise a total of 507 006 datapoints covering 
50 complete transits. 58 161 of these datapoints are close to tran- 
sits, and I phase-binned them by a factor of 50 to get 1164 nor- 
mal points. I ado pted e cos uo = 0.000 6 ± 0.0099 and esinoj = 
0.0350±0.0170 iBuchhave etal .1201 ih . whch suggests that eccen- 
tricity is significant at the 2a level. The RP errorbars were smaller 
than the MC ones, and I was able to adopt LD-fitted solutions (Ta- 
ble A52). 

Given the extremely high quality of the light curve I tried fit- 
ting for 1/3 directly. Investigations of this possibility in Paper III 
showed that third light is almost completely degenerate with other 



photometric parameters (rA, and i), so may not be obtainable 
even from exceptionally good data. I obtained = 0.223 ± 0.085 
which is within Icr of the expected value but still poorer than us- 
ing outside constraints. I also checked whether it was po ssible to 
fit for ecosw and/or esinoj, and confirmed the result of iKippin^ 
(2008') that such parameters cannot be extracted. The best-fitting 
values were ecoscj = 0.350 ± 0.065 and esintj = 0.149 ± 
0. 019, which are inconsi stent with the known results determined 
by IBuchhave et alj ( 1201 ih from RVs and a Spitzer occultation. 

A summary of the light curve solutions is given in Table A53 
and the best fit is plotted in Fig.[T2] I find somewh at diff erent pho- 
tometric parameters from those of Buc hhave et al] ^201 ih . which is 
likely down to the use of short-cadence data in the current work. 
The third light has a strong effect on the results, as the photometric 
parameters are all strongly correlated with it. It thus sets a limit on 
the quality of the solution, so improved measurements of Ls would 
be the best way to improve our understanding of the system. 

I find a rather different set of physical properties for Kepler- 
14 (Table A54) , with a le ss massive but more evolved star. Com- 
pared to Buchh ave et al.l (|2011) I find a stellar mass of A/a = 
1.318 ± 0.052 ± 0.029 M© rather than 1.512 ± 0.043 Mq (a 
change of 2.6a), a very similar stellar radius and a 1.4o" dif- 
ferent surface gravity (logpA = 3.918 ± 0.040 ± 0.003 ver- 
sus 3.994ln n3 fi). The log pa deri ved from spectral analysis is 
4.11 zb 0.10 ( Buchhave et a mm J). These changes propagate to 
the planet, which is found to be less massive but the same size as 
before. Kepler will provide further photometry to help investigate 
this system, but obtaining revised T^a, [§] and L3 values would 
be the best way to improving our understanding of the Kepler-14 
system. 
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5.13 Kepler-15 

lEndl et al] ( 1201 ih presented the discovery of this planetary system 
containing a comparatively low-mass (0 .70 ± 0.10 Mjup) planet 
around a cool (5595 ± 120 K) star. Endl et al.l r201 i h used the Ke- 
pler long-cadence photometry from quarters Ql to Q6 plus RVs 
from two telescopes to characterise the system. However, Kepler 
short-cadence data are also available for quarters Q3 to Q6, and 
their much greater time resolution should yield better results than 
for the long-cadence data. 

The Kepler Q3 to Q6 short-cadence data were modelled in- 
dividually quarter-by-quarter, using the third light values from the 
KIC and adopting an uncertainty of 0.5% in these va l ues. E ccen- 
tricity was not found to be significant by lEndl et al.l ( 1201 ih so I 
assumed e = 0. Each dataset was phase-binned by a factor of 20, 
giving 613, 563, 685 and 644 normal points respectively. The full 
solutions are in Tables A55 to A58 and the best fits are shown in 

Table A59 summarises the light curve solutions, and shows 
that their interagreement is good (except for k where = 1-8)- 
The final photometric parameters a re extremely prec ise. However, 
their ageement with the findings o fl Endletal.ll201lh IS very poor: 
differences of 24cr occur for ta and i, 7a for k and 33(7 for rb. 
It is unclear how such a large discrepancy could arise, although 
one possibility to investigate would be the use of long-cadence 
ve rsus short-cadenc e data. I was unable to reproduce the results 
of lEndletalJ ^1 l|) by fitting the long-cadence data either with 
or without numerical integration. My new results should be pre- 
ferred because the short-cadence data are much better for measur- 
ing planet parameters than the long-cadence data (see Paper IV). 

As would be expected from the rather different photometric 
parameters, my calculated physical properties (Table A60) are at 
odds with those of Endl et al. (2011). The most obvious changes 
are for the radii of the two components, wh ich go from Ra = 
0.992ta"om R© and i?b = 0.96t^;^? Rjup fendl et al.ll20"llh to 
Ra = 1.253 ± 0.047 ± 0.020 R© and i?b = 1.289 ± 0.050 ± 
0.021 Rjup (this work). The loggA I find is in better agreement 
with the spectroscopic measurement, albeit at low significance. My 
revisions to the physical properties of the Kepler-15 system de- 
crease the planet's measured density by more than a factor of 2, 
increase the measured age of the star, and include the first mea- 
surement of gb, Tcq and O. 



5.14 Kepler-17 

Kepler- 17 comprises a hot and massive planet (Tj^ — 1712 ±25 K, 
Mb = 2.34t«:0| Mjup) orbiting an active star. Starspot activity is 
obvious both in the overall light curve and by spot crossings during 
the planetary transits, allowing a measurement of t he star's rota- 
tion period in the discovery paper toesertetal.l201ll) . The rotation 
period is close to eight times Porb, giving a 'stroboscopic' effect 
whereby some starspots can be identified each eighth transit. This 
also allows the orbital obliquity to be constrained to be less than 
15° toesert et al.ll201 ll) . Spitzer observations of an occultation al- 
lowed limits to be placed on orbital shape e < 0.011 in combina- 
tion with the Kepler photometry and RV observations. An indepen- 
dent d iscovery of the Kepler-17 system was ma de bv.Bonomo et al.. 
( |2012|) . who combined their RVs with those of lPesert et alj j201lh 
to obtain a more precise Ka- Their other result were in good 
agreement with those from the discovery paper, but less precise 
as they only had K epler Q1-Q2 data versus the Q1-Q6 available to 
lDesertetal.l ( l201lh . 
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Figure 14. The Kepler Q4-Q6 short-cadence light curves of Kepler-17 com- 
pared to the JKTEBOP best fits. Other comments are the same as in Fig.|6] 

As before (see HAT-P-1 1 in Paper IV) I have assumed that the 
starspot activity under the transit chord is on average the same as 
that for the rest of the stellar photosphere. This means that the data 
affected by starspots do not require special treatment, and the effect 
of the starspot crossings on the solution will simply be to inflate the 
errorbars as an additional source of correlated noise. This in turn is 
partially coverted into white noise when the data are phase-binned, 
and will ultimately cause larger errorbars to be found than if the 
star were quiet. 

The Kepler Q4 to Q6 data are available in short-cadence 
flavour, and these were phase-binned as usual to yield 621, 634 
and 621 normal points, respectively. A significant fraction (~15%) 
of transits were normalised using a quadratic function instead of a 
straight line, in order to cope with the modulation due to starspot 
activity. 120 datapoints were rejected from the Q5 data by a 4(t 
clip, due to an unusually large number of scattered observations, 
whereas the numbers rejected from Q4 and Q6 data were 3 and 1 1 
respectively. 

Third light was tackled by adopting the season-dependent con- 
tamination values from the KIC and assigning an uncertainty of 
0.5% to each one. The orbit was assumed to be circular. The full 
sets of light curve solutions are given in Tables A61 to A63 and in 
each case the RP errors were found to be moderately larger than 
the MC ones, due to the starspots on the stellar surface. The best 
fits are plotted in Fig.[T4l where the influence of spot activity can 
be seen in the residuals of the best fit. The light curve solutions are 
summarised in Table A64, and are in good agreement, except for 
k where Xu = 5.8. This is due both to starspot activity and the 
need to use quadratic functions to normalise some transits, and has 
been accounted for in the errorbars. It is not sufficient to cause a 
disagreements in ta or rb. 

The light curve solution I find is close to a central transit 
deg) whereas [Pesert et al] ( 1201 l l) prefer a lower- 
inclination solution {i = 87.2 ± 0.15 deg) which differs by 2.6(t 
from my own. This affects ta (which is 3.9(t larger according to 
(Desert et al.i) and k (3.9a" smaller). The source of this divergence is 



© 0000 RAS, MNRAS 000, 000-000 



14 John Southworth 



1.000 



0.995 - 



S 0.990 



0.985 



0.980 




-0.04 



-0.02 



0.00 
Orbital phase 



0.02 



0.04 



Figure 15. The Kepler light curves of KOI-135 compared to the JKTEBOP 
best fits. The upper dataset is the Q1-Q3 long cadence data and the lower 
dataset is the Q3 short cadence data. Other comments are the same as Fig.[T] 



unclear. iBonomo et alj ( I20l3) had access to only a litt le of the Ke- 
pler d ata, so adopted the photometric parameters from lOesert et al.l 
( |20T l|) rather than calculate their own. Given the discrepancies, the 
larger errorbars found in the present study should be preferred to 
those given in the literature. 

When calculating the physical properties of the Kepler- 17 sys- 
tem I found slight disagreements within results using the five sets of 
stellar models, which leads to a significant systematic component in 
the errorbars of the final system parameters. The values themselves 
are in reasonable agreement with previous studies (Table A65). As 
Kepler continues to gather photometry of this system, the action 
with the most benefit would be the procurement of improved spec- 
troscopy for Tea , and RV measurements. 

5.15 KOI-135 

KOI-135 was discovered bv lBonqmo et al. ( 20121) to co ntain a mas- 
sive planet (3.09 ± 0.21 Mjup). Bonomo et al. j2012h had access 
to only the Ql and Q2 long-cadence data. One of the problems 
they faced is that the Kepler long-cadence sampling rate is close to 
an integer multiple of the Porb of KOI-135, leading to very poor 
sampling of the shape of the transit (see Fig.|15b. Q3 short-cadence 
data have now been released, which should allow the current work 
to produce a significant improvement in the measured system pa- 

ram eters. 

iBonomo et al.l ( 1201 2h found that their spectroscopic measure- 
ment of log (?A disagreed with that obtained from analysis of the 
transits by 2.7cr. They found that using the transit-derived loggA 
resulted in an insignificantly different measurement of T^a from 
their spectra. 

My first action was to obtain an improved orbital ephemeris 
using the Ql to Q3 data. This was done in the same way as for 
Kepler-7 in Sect. l5.10l resulting in: 

To = BJD(TDB) 2 455 068.235237(54) + 3.0240933(27) x E 

The precision of this Pp rb is a factor of ten b etter than that based 
on only Ql and Q2 data (Bonomo e t al.ll2012h . 

Constraints on third light were obtained from the KJC as 
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Figure 16. The Kepler light curves of KOI-196 compared to the JKTEBOP 
best fits. From top to bottom the three datasets come from the Ql, Q2 and 
Q3 time period. Other comments are the same as Fig.[T] 



above, and errors of 0.5 % specified. A circular orbit was adopted 
as Bonom o et al.l ( 1201 2h found e < 0.025 (significance level not 
stated). Only the short-cadence (Q3) data were modelled, for the 
reasons given above. The original 125 140 datapoints were slimmed 
down to 18 938 close to any of the 27 transits observed by Kepler, 
13 were rejected by a Aa clip, and the remainder were phase-binned 
by a factor of 25 to get 759 normal points. Table A66 shows that the 
LD-fitted solutions can be used. The RP errors were smaller than 
the MC errors. The Q3 d ata yield a Scr smaller r a and rb than found 
by IBonomo et alj ( l20I2h . so the availability of improved photome- 
try has resulted in the modification of the measured system proper- 
ties (Table A67). 

Despite the large uncertainties in Tcff (143 K) and [^] 
(0. 1 1 dex), I found highly consistent physical properties for the 
KOI-135 system using different sets of stellar models (Table A68). 
Both the star and planet are slightly smaller than found by 
IBonomo et al.l ( l2012h . but these changes are within the errorbars. 
Forthcoming Kepler photometry and further spectroscopic mea- 
surements of Tcff and [^] would be very helpful in refining the 
system properties. 



5.16 KOI-196 



ISanteme et al.l (1201 lal) discovered the planetary nature of this sys- 
tem, based on Q1-Q2 data f rom the Kepler satellite and RVs from 
the SOPHIE spectrograph. ISanterne et al.l also detected occulta- 
tions in the Kepler data, of a depth which suggest the planet has 
a high albedo and at an orbital phase which is consistent with a cir- 
cular orbit. Since this study the Q3 data have become available, but 
like Ql and Q2 are only long-cadence. 

Third light was estimated from the KIC, 0.5% errors assigned, 
and a circular orbit assumed. An improved orbital ephemeris was 
obtained from the Q1-Q3 data using the method in Sect. |5.10l 

To = BJD(TDB) 2 455 066.669310(48) + 1.8555588(15) x E 
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Figure 17. The Kepler Ql, Q2 and Q3 long-cadence light curves of KOI- 
204 compared to the JKTEBOP best fits. Other comments are the same as 
Fig.[I] 

For the Ql data, 257 of the 1625 datapoints are in the region 
of one of the 18 transits. For Q2 these numbers are 4802, 648 and 
45, and for Q3 they are 4145, 662 and 46. The data were not binned 
or phased, and TVint = 10 was used in the modelling process. Two 
families of solutions were found: one with a nearly central tran- 
sit (i Ri 90°) and one with i « 80-82°. The i ^ 80-82° family 
occurs mainly for the LD-fixed light curve solutions, and results 
in weird physical properties. In contrast, the LD-fit/fix solutions 
universally have i « 90° and lead to reasonable physical proper- 
ties. This implies that fixing LD coefficients can have a deleteri- 
ous elfect on a light curve solution, specifically when working with 
long-cadence data. I adopt ed the LD-fit/fix solut ions (Tables A69 
to A71), in agreement with lSanteme et al] ( 1201 lah . and found that 
the RP errors are of similar size to the MC errors. 

The photometric parameters are given in Table A72 and the 
best fits are plotted in Fig.[T6] The interagreement between the 
Q1-Q3 data is excellent, but further photometry (pr eferably short- 
cadence ) is needed to obtain a definitive solution. iSanterne et al.l 
( l2011al) found similar but slightly discrepant photometric parame- 
ters, based on fewer data. 

As with KOI- 135, I found highly consistent physical prop- 
erties for the KOI- 135 system from different sets of stellar mod- 
els (Table A73) despite the lar ge uncertainties in Teff (100 K) and 
[^] (0.11 dex). Compared to ISanteme et al.l j2011a l). I deduced 
a slightly smaller and less massive star and planet. Further spec- 
troscopy and forthcoming Kepler photometry would allow a useful 
improvement in measurements of the physical properties of KOI- 
196. 

5.17 KOI-204 

KOI-204 b was discovered by iBonomo et al.l | |2012|) as part of a 
program t o observe the brighter Kepler TEP candidate systems with 
SOPHIE. iBonomo et al.l ( l2012h had access to Kepler Ql and Q2 
data, and Q3 data are now available. I derived an improved orbital 




-0.03 -0.02 -0.01 0.00 0.01 0.02 0.03 
Orbital phase 



Figure 18. The Kepler Ql to Q5 light curves (from top to bottom) of KOI- 
254 compared to the JKTEBOP best fits. Other comments are the same as 
Fig-E] 

ephemeris using these data and the method in Sect. |5.1()] finding 
To = BJD(TDB) 2 455 067.02680(15) + 3.2467220(73) x E 

The orbit was assumed to be circular dSonomo et al.l found 
e < 0.021) and third light values were obtained from the KIC. 
These are a little larger than usual, being 0.060, 0.093, 0.154 and 
0.044, in seasons to 3 respectively. The individual quarters of data 
were then fitted using TVint = 10, and the full solutions are given 
in Tables A74 to A76. The RP errors are smaller than those from 
the MC algorithm. The best fits are plotted in Fig.|17| 

Table A77 gives the final photometric parameters. The three 
datasets agree well except for k {xi ~ 1-9), which is likely a symp- 
tom of spot crossings duing the transit. These are not visible to the 
eye, due to the poor time resolut ion of the long-cadence data. The 
values from lBonomo et al] ( 1201 2h are concordant with my own, but 
their errorbars are smaller despite being based on significantly less 
data. The same comments apply to the physical properties deter- 
mined and gathered in Table A78. More observations of all types 
are needed to characterise this planetary system well. 

5.18 KOI-254 

Ijohnson et al. discovered this planetary syst em from RV 

follow-up of Kepler KOI objects dBorucki et alj201ll) classified as 
M dwarfs, in order to probe this host star mass regime. The stellar 
component is on the border of the M spectral class (Ma = 0.57 ± 
0.06 Mq and T^s = 3820 ± 90 K), and the planet is the largest 
one found around a sub-4000 K star (7?b = 1.00 ± 0.07 Rjup) by 
some distance (the next is GJ 436 b at 0.37 ± 0.02 Rjup). KOI-254 
is therefore an interesting object in delineating the minimum stellar 
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mass required to host a Jupiter-sized object. I Johnson et alj ( 1201 2h 
analysed the Kepler Ql and Q2 data, which shows clear spot activ- 
ity, plus a Z-band transit light curve obtained using the 1 m Nickel 
telescope at Lick Observatory. 

At the time of writing, Kepler Ql to Q5 data were available, 
all at long cadence. A refined orbital ephemeris was calculated as 
in Sect. lSlOl 

To = BJD(TDB) 2 455160.956396(31) + 2.45524122(65) x S 

which is a factor of ten more precise than the ephemeris based on 
only the Ql and Q2 data. 

The data from each quarter were modelled assu ming an eccen- 
tric o rbit with e = O.lllo;^ and u = 230° ± 68° ijohnson et all 
l2012t) . third light values from the KIC with uncertainties assumed 
to be 0.5%, and iVint = 10 (Tables A79 to A83). Table A84 shows 
the agreement to be very good both betw een the five datasets and 
with the values from I Johnson et al] ( |2012|) . The RP errors were all 
smaller than the MC errors. The best fits are shown in Fig.1181 

Precise physical properties become difficult to obtain into the 
M-dwarf regime, as the interagreement between stellar models de- 
teriorates, but KOI-254 is thankfully not strongly affected. The 
physical properties I find (Table A85) agree well with those of 
Ijohnson et al.l ( |2012|) . but the availability of additional data allows 
the errorbars to be beaten down. Kepler will deliver further pho- 
tometry, but more spectroscopy would be useful in nailing down 
the Tefi, [^] and spectrosocpic orbit of the host star. A significant 
source of uncertainty in the current results stems from the mea- 
surement of e and uj, which is rather imprecise due to the small RV 
signal exhibited by the star. 

5.19 KOI-423 

KOI-423 contains an object of brown-dwarf mass (Mb = 17.9 ± 
1.8) in a long-period (21.1 d) eccentric orbit (e = 0.12llo;{523) 
around a hot and m etal-poor (Tpg = 6 260 ± 140 K and [^] = 
-0.29 ± 0.10) star ( iBouchv et alj201lh . Its mass, however, places 
it firmly in the planetary ra ther than the stellar population according 
to lOrether & Lineweaver (2006). Bouchv et al. (2011) based their 
analysis on RVs from the SOPHIE spectrograph, and the Kepler Ql 
and Q2 data (which encompass only five transits between them). 

The Q3 to Q6 data are also now public; all six datasets are 
not only long-cadence but also sample the transit shape poorly 
(Fig.|19t. A revised orbital ephemeris was first procured using the 
method in Sect. l5.10l 

To = BJD(TDB) 2 455 204.55515(14) + 21.087168(21) x E 

Each dataset was solved with numerical integration (A'^int = 
10), constraints on L3 from the KIC, and adopting e cos lj = 
0.120l[J:g24 and esino; = -0.019lE!Ui| dSouchv et al]|201lh . 
A quadratic or cubic function was needed to normalise the ma- 
jority of the transits to unit flux. MC errors were larger than RP 
errors, and full solutions (except for LD-fitted) can be found in Ta- 
bles A86 to A91. The total dataset covers 22 transits, of which one 
had to be rejected due to normalisation problems caused by miss- 
in g data, which is a si gnificant improvement on the five available 
to lBouchvetai]j201lh . 

Table A92 shows that the six different light curves agree well. 
The combined solution has a somewha t larger i, va and V h, but 
these are within l-2a of those found bv lBouchv et al.l ilOl ih . The 
physical properties of KOI-423 (Table A93) are tole rably well de- 
fined, and are reasonably close to the values found bv lBouchv et al.1 
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Figure 19. The Kepler Ql (top) to Q6 (bottom) light curves of KOI-423 
compared to the JKTEBOP best fits. Other comments are the same as Fig.[T] 

( I2OI ll) . Kepler short-cadence data would be valuable for this object, 
as would further spectroscopy to refine Tag and [|^] . 

5.20 KOI-428 

This planetary system was discovered by ' Santeme et al. I( l2011bh . 
and studied using the Kepler Ql data. It was subsequently analysed 
in Paper IV, based on the Ql and Q2 data. I have revisited it because 
Q3 to Q6 data are now available - all long-cadence - and cover 67 
transits compared to the four in Ql and the eight in Q2. The other 
modification to my previous work concerns the small L3, which in 
Paper IV was ignored but in the current work is taken from the KIC 
as usual. The new orbital ephemeris based on the Ql to Q6 data is: 

To = BJD(TDB) 2 455 204.83996(20) + 6.8731697(96) x E 

A circular orbit was assumed, and each quarter of data was 
tackled separately with Mnt = 10. The full solutions are given in 
Tables A94 to A99 and the best fits displayed in Fig.[20] Correlated 
noise is unimportant, and the six light curves agree within the rather 
large MC errorbars (Table AlOO). This underlines the low informa- 
tion content in data with a sampling rate of only one per 30min. 
Short-cadence data is necessary to obtain good results on this ob- 
ject. The revised ta is comfortably larger than found in the two 
previous studies; the change is within the errors for Paper IV bu t 
outside the rather small errorbars quoted bv lSanterne et alj ( l2011bl) . 

The revised physical properties (Table AlOl) are significantly 
different, due to their sensitivity to ta via the stellar density. 
The masses are almost unchanged but the radii both increase by 
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Figure 20. The Kepler Ql (top) to Q6 (bottom) light curves of KOI-428 
compared to the JKTEBOP best fits. Other comments are the same as Fig.[T] 



about 10%. The star is now quite evolved (loggA ~ 3.812 ± 
0.048 ± 0.017) and the largest known TEP host star at Ra = 
2.48 ± 0.17 ± 0.20 R© (the second-largest being Kepler-14A at 
2.09 ± 0.11 ± 0.02 R©; Sect. l5TT2] l. A significant contribution to 
the uncertainties in these numbers comes from significant disagree- 
ment between the five model sets. The DSEP and 7^ models pre- 
fer a larger system scale and a slightly pre-main-sequence star, 
whereas the other three models find an age of 1-2 Gyr. Arbitrarily 
restricting the star to be on the main sequence would improve the 
agreement between models, but cannot be scientifically justified. 
Further Kepler photometry and ground-based spectroscopy would 
be important steps towards confirming and refining the physical 
properties of this planetary system. 



5.21 OGLE-TR-56 

OGLE-TR-56 has been studied in previous papers of this series, but 
qualifies for further attention due to the recent publication of much 
improved photometric d ata. Paper I summaris es the observational 
history of OGLE-TR-56. lAdams et al] ( 1201 Ibl) have since obtained 
19 new transit light curves, primarily with the high-resolution 
Magic camera on the Magellan telescopes. A few transits were also 
observed using the POETS and IMACS ima gers on the same tele - 
scopes, and with GMOS on Gemini-South. lAdams et al. I( l2011bl) 
fo und a discrepancy between their transit duration and that obtained 
bv lPont et alj ( |2007|) from one transit observed using VLT/F0RS2. 
A re-analysis of the VLT data resulted in a smalle r discrepancy 
but more scattered photometry. One observation bv lAdams et al.l 
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Figure 21. The lAdams et alj jlOllbh hght curves of OGLE-TR-56 com- 
pared to the JKTEBOP best fits. The sources of the data and the filters used 
are annotated on the left and right, respectively. Other comments are the 
same as Fig.|6] 



is of the same transit observed using the VLT, and the 
transit shape is in g ood agreement with the others presented by 
lAdams et alj ( 1201 Ibl) . I have therefore discounted the VLT data (as 
analysed in Paper I) as possibly suffering from systematic noise, 
and rederived the pr operties of OGLE-TR -56 based on the pho- 
tometry obtained bv lAdams et ai]( l2011bh . The datasets obtained 
using the SITe CCD and B, r and z filters were ignored as each do 
not cover a whole transit. 



5.21.1 Light curve analysis 

The data taken bv lAdams et"al] ( 1201 Ibl) using the SITe CCD and 
/ filter cover four transits, of which one was only partially ob- 
served and suffers from correlated noise so was rejected. 498 data- 
points cover the remaining three transits and were fitted with Porb 
fixed. The RP errors are slightly larger than the MC errors, and the 
LD-fixed solutions were adopted because the LD-fit/fix alternatives 
gave unphysical LD coefficients (Table A102). 

Two transits were observed using POETS in the Is band (a 
filter described as 'Schuler Astrodon lohnson-Cousins'). There is 
a short-lived increase in flux just before second contact for the first 
transit, which could be due to a stellar flare. The observations were 
taken at a high cadence, resulting in 8961 datapoints. These were 
binned in time by a factor of 20 to obtain 449 normal points, with 
four datapoints rejected by a 3a clip. Once again, the RP errors 
were slightly larger than the MC ones, and the LD-fixed solution 
had to be adopted (Table A103). 
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One more transit was observed using POETS without a fil- 
ter, in an attempt to achieve a lower P oisson noise. Thi s tran- 
sit event was the same one observed by IPont et aD ( l2007h using 
VLT/F0RS2. The observations were taken at high cadence, result- 
ing in 3245 datapoints. Those after time HJD 2453936.68 were re- 
jected as they are well after the completion of the transit and suffer 
from systematic noise. The remainder were time-binned by a fac- 
tor of ten to give 95 normal points, then modelled using JKTEBOP. 
LD coefficients appropriate for the Cousins R band were used. The 
LD-fit/fix solutions are unfortunately poor (Table A104), and the 
LD-fixed solutions are possibly biased due to the imposition of R- 
band LD coefficients onto an unfiltered dataset. I therefore did not 
consider these data when calculating the final photometric parame- 
ters. 

One transit (352 datapoints) was observed using IMACS and 
a wide-band filter (630-950 nm), for which I adopted Cousins I- 
band LD coefficients. Correlated noise was unimportant, but once 
again LD-fixed solutions had to be adopted (Table A105). Another 
transit (307 datapoints) was obtained using Gemini/GMOS and an 
i filter. For this one correlated noise was unimportant and it was 
possible to adopt the LD-fit/fix solution (Table A106). 

Finally, nine complete transits were observed in the i band us- 
ing Magellan/Magic and its E2V CCD. For these 5186 datapoints 
I found that correlated noise was moderately important but that it 
was also possible to adopt the LD-fit/fix solution (Table A107). All 
light curves are compared in Fig.|2T]to their respective best-fitting 
models. 

One minor concern is that OGLE-TR-56b is sufficiently close 
to its parent star to be distorted beyond the formal limit of applica- 
bility of the EBOP model. All other work on OGLE-TR-56 has been 
pe rformed using spherical geo metry, such as th e analytical models 
bv lMandel & Agoll ( [20o3) and lGimenel ( |2006|) . JKTEBOP is more 
sophisiticated than those, in its use of biaxial spheroids to repre- 
sent the stellar and planetary figures, but it would be worthwhile to 
check them using a more physically realistic model which incorpo- 
ra tes Roche geo metry. An excellent discussion of this can be found 
in lBudail ( I201TI) . 

The agreement between the six different datasets from 
I Adams et al. (2011b) is very good, ranging from Xv ~ 0.19 for i 
to Xv = 0.98 for k. The combined solution was calculated exclud- 
ing the unfiltered POETS data, and the errorbars were (as usual) 
conservatively not reduced in light of the low Xv values. The solu- 
tions are summarised in Table A108, which also hosts a comparison 
with (the many) previous results for OGLE-TR-56. M y new photo- 
metric p arameters are in good agreement with those of lAdams et al.l 



metric p ; 
( l2011b[) 



but are in severe disagreement with all previous work. This 
is expected because my new results are based on the data from 
lAdams et alJboilbl . whereas all prior analyses had access to only 
the original survey-quality OGLE data and/or the VLT observations 
which have turned out to be affected by red noise jAdams et al.l 
l2011bl) . 



5.21.2 Physical propertie.s 

A major change in the photometric parameters induces a similar 
effect in the physical properties of the system. My results (Ta- 
ble A 109) show that there is some disagreement between differ- 
ent model sets, in that the Claret and Teramo models coverge on 
solutions with a la lower Tea. This causes a significant system- 
atic uncertainty to enter the error budget, but does not obscure the 
main story. I have found physical properties which are very differ- 
ent from all previous estimates, and point to OGLE-TR-56 A be- 
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Figure 22. Mass-radius plot showing the known TEPs (blue squares for 
those in the Homogeneous Studies series and red triangles for the remain- 
der) with their errorbars suppressed for clarity. The new result for OGLE- 
TR-56 b is shown with a large black filled circle, whereas previous results 
are shown with small black filled circles. 



ing somewhat evolved (Ma = 1.339 ± 0.077 ± 0.066 and 
i?A = 1.737 ± 0.035 ± 0.029 R0). By contrast, all nine pre- 
viously published measurements of Ra were in the range 1.10- 
1.36 Rq. The same situation is found for the planet (Mb = 
1.41±0.17±0.05 Mjup), for which the new radius measurement is 
Rh — 1. 734 ± 0.051 ± 0.029 R,7up whereas all previous values are 
in the range 1.23-1.378 R,Tup. This is illustrated in Fig.[22] OGLE- 
TR-56 b therefore joins the group of the most inflated TEPs. This 
is unsurprising insofar as its host star is relatively hot and mas- 
sive. A detailed discussio n of the most inflated TEPs is given by 
ISouthworth et alj ( l2012bh . 

How did such a large chang e in the measured properties oc- 
cur? Early analyses of the svstem ([Konacki et alj2003l ; [Torres et"ail 
|2004 ISantos et al.ll2006l : iBouchv et al.ll2005h all relied on on the 
OGLE discovery light curve ( IUdalskietalj2002bl) which is of low 
qualit y compared to d edicated foll ow-up observations. Later anal- 
yses ( Pont et alj2007l : [Torres et a l.' 2008: Paper II; Paper III) rested 
on the V LT light curve from IPont et al.. t2007i) . which has been 
shown by lAdams et alj 1 I2OI Ibh to suffer fro m substantial system- 
atic noise. It is only with the data gathered bv I Adams et al. Il l2011bl) 
that it is now possible to derive definitive physical properties of the 
system. 

lAdamsetal.U20Ilbl) themselves studied their light curves in 
detail, but adopted old values of Ma and Ra from Paper III for cal- 
culating the system properties. They therefore did not propagate the 
stellar density found from their light curves to the physical proper- 
ties of the system, and so ended up with properties similar those 
found in previous works. It is only when the newer photometric 
data are used to specify the stellar density, and therefore its mass 
and radius in conjunction with theoretical models, that the system 
parameters experience a significant change. The physical proper- 
ties of OGLE-TR-56 found in the current work are the first ones to 
provide a full picture of the system. 

One moral from this story is that a single good light curve 
is not enough to provide a definitive set of physical properties, as 
undetected sysmatic errors may bias the results. A significant frac- 
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Figure 23. The Magell an jAdams etai]l2010ah and VLT jMinniti et all 
I2OO7I: iHover et al.ll201 ll) light curves of OGLE-TR-111 compared to the 
JKTEBOP best fits. Other comments are the same as Fig.|6] 



tion of oth er TEPs might be sub ject to this con cern, for example 
HP 8 0606 ( lHebrardetai]|2010l) and WASP-7 jSouthworth etal] 
I2OIII) . A second moral is that any significant revision to the photo- 
metric parameters should be propagated into the physical properties 
of the system. 



5.22 OGLE-TR-111 

The transitin g natur e of OGLE-TR-111 was announced by 
Udalski et alj j[2002ah and its planetary nature confirmed by 



Pont et alT i 2OO4I) . Its publication hi story was summa rised in Pa- 
perl (sect. 3.9), where the data from lWinn et al.l ^20071) were anal- 
ysed. Three new sets of transit light curves have be en published 
since th en, prompting this rediscussion of the system. [Adams et al] 
( 12010 a') observed six transits through an SDSS i filter, using 
Magellan/MagIC-e2v high-resolution imager, allowing them to im- 
prove the measurement of Rh and d i smiss t he existence of the 
TTVs tentatively found by iDfaz et al.l j2008l) . iHover et al.l j201lh 
studied five transits, one with VLT/FORSl and a "v-HIGH" filter 
and the remaining four with VLT/FORS2 and a Bessell / filt er. I 
have a lso included the single l/-band light curve from^inniti e t al.l 

( l2007h in my analysis. 

The six transits observed by I Adams et al. I ( l2010ah are con- 
tained in 2205 datapoints, for which I found that correlated noise 
was moderately important. I fitted for Porb, and the li ght curve 
solutions are given in Table Al 10. iMinniti et al] ( l2007h obtained 
213 datapoints covering one transit in V, for which I f ound no evi- 
dence of correlated noise (Table Al 1 1). The data from lHover et al.l 
( I20T ij) were supplied with timestamps in JD, which I converted to 
BJD(TDB) before analysis. The V-band transit has 323 datapoints 
and the four J-band transits comprise 1967 datapoints. I found the 
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Figure 24. The Magellan light curve of OGLE-TR-113 <Adams et al] 
IjOlOb) compared to the JKTEBOP best fit. Other comments are the same 
as Fig.[6] 



RP errors to be larger than the MC errors for both datasets (Tables 
A112and A113). 

The light curve solutions are summarised in Table Al 14 and 
the best fits are shown in Fig.|23] The agreement between the light 
curves is good, except for k which has ~ 1-9- This situation 
has occured frequently within this series of papers. Published pho- 
tometric properties agree well with my own. A similar story occurs 
for the physical properties (Table Al 15), although I find a com- 
paratively low concordance between the results for different stellar 
models. OGLE-TR-1 1 1 would benefit from an improved measure- 
ment of A'a, although its faintness {V — 16.96) makes this non- 
trivial. 



system by 



5.23 OGLE-TR-113 

OGLE-TR-1 13 was also found to be a transitin g 

lUdalski et in i2002ah . and a planetary system by iBouchv et al.l 
(L2OO4I). It was studied in P aper III (sect . 4.3 and fig. 7), base d on the 
photometr y presented bv [Gillon et al . ( 2006), Snellen & Covind 
(I2OO7I) and lDiaz et al] ( |2007|) . Since"then l Adams et alj ( 1201 Obb have 
published six nice transit light curves from Magellan/MagIC, al- 
lowing them to revise J?b and set limits on TTVs. They found a 
possible gradual decrease in Porb, with a significance of 4(t, which 
should be investigated onc e a longer temporal b aseline is available. 

I fitted the data from lAdams et alj d2010bl) . which represents 
2697 datapoints, including Porb as a fitted parameter. My initial 
model was compromised by some poor data which w e re inclu ded 
in the datafile but not in the analysis by Adams et al. (2010b). 79 
datapoints were rejected, allowing a good fit to be found to the 
remaining data. I found correlated noise to be slightly important. 
The LD-fit/fix solution w as adopted (Table A 1 1 6 and Fig.l24b. 

The new data from lAdams et al. I ( l2010hl) allow a great im- 
provement in the photometric parameters (Table Al 17) and in the 
physical properties (Table All 8). I found a good agreement be- 
tween results from different model sets, and against values from 
the literature. As with OGLE-TR- 1 1 1 , a more precise if a would 
be useful but difficult to achieve because the system is faint (V = 
16.48). 



5.24 OGLE-TR-132 

lUdalski et alj j2003h observed transits in this system, and 
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Figure 25. The Magellan light curves of OGLE-TR-132 from lAdams et alJ 
i2011al) . compared to the JKTEBOP best fit. Other comments are the same 
as Fig.|6] 



iBouchv et alj ( l2004h confirmed its planetary nature. OGLE-TR- 
132 was studies in Paper I (s ect. 3.10 and fig. 12 ) using a nice 



VLT/FORS2 light cur ve from iGillon et al.l ( l2007h . I have revis 



ited the system because IXdams et alj ( 1201 lah has published a light 
curve of 4940 datapoints covering seven transits, using them to 
measure a revised Rh and put upper limits on the presence of TTVs. 
Two of these transits were observed at a much higher cadence (10- 
15 s) versus the rest (35-130 s), so these were time-binned by a 
factor of five. This is an important step if the data are modelled 
together and subjected to an RP algorithm. 

The resulting 2131 datapoints were analysed, with the find- 
ings that modest correlated noise is present and that the LD-fit/fix 
solution is good (Table Al 19 and Fig.|25l). Final photometric pa- 
rameters were calculated from these results and those in Paper I. 
The agreement between these results and also with published val- 
ues is excellent (Table A120). The ensuing physical properties (Ta- 
ble A121) are an improvement over Paper I due to the additional 
data, but are in good agreement with literature results. OGLE-TR- 
132 needs more RVs if its measured physical properties are going 
to be improved significantly. 
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Figur e 26. The GROND light curves of OGLE-TR-L9 from iLendl et alj 

1 2010) compared to the JKTEBOP best fits. Other comments are the same as 
Fig.! 



from lSnellen et al.l ( l2009h or Paper III (both of which used the same 

sma ller dataset). 

iLendl etall feOld) did not derive the full physical properties 
of OGLE-TR-L9, so Table A127 contains the first physical proper- 
ties measured using their extensive photometry. The parameters are 
consistent with previous solutions, but with a useful lowering of the 
errorbars due to the new light curves. An improved Ka and [^] 
are needed; the light curves are unusually not the observational bot- 
tleneck because these spectroscopic constraints are so poor (mainly 
due to the faintness of the system at F « 15.5). 



5.25 OGLE-TR-L9 

This was discovered by ISnellen et alj j2009l) based on archival 
OGLE data, follow-up four-band photometry from the GROND 
imager, and UVES spectroscopy. It is of interest as containing 
a massive planet (4.4 ± 1.5Mjup) around an unusually hot star 
(6933±58 K), and was duly studied in Paper III. lLendl et aljfioTol) 
have recently obtained five more transits with GR ONEj I have 
model led t he combined datas et of six transits from ISnellen et al] 
( l2009l) and lLendl et al. 112010') as they were obtained with the same 
instrument and observing strategy. The numbers of datapoints is 
as follows: 922 in g, 921 in r, 917 in i and 922 in z. In all cases 
the RP errors were smaller than the MC errors. The full sets of 
solutions can be found in Tables A122 to A125, the summary in 
Table A126, and the best fits in Fig.|26l The final photome tric pa- 
rameters agree well with those found bv lLendl et alj ( l2010i) . as ex- 
pected given the identical datasets used, but not so well with those 



^ The datafiles as lodged on CDS are incorrect but the data are available on 
request to M. Lendl. 



5.26 TrES-4 

This system, containing a large and low-density planet, was stud- 
ied in Paper III based on the ph otometry presented in the discovery 
paper (Mandush ev et alj|2007l) . New data are now available from 
Chanetal. (20lC covering five transits in the i band obtained us- 
ing KeplerCam. Here I present revised results incorporating these 
new data. A notable result since Paper III is the me asurement of 
the Rossiter-McLaughlin effect bv lNarita et ai] l l2010l) . showing the 
system to be aligned to within l.Scr. 

TrES-4 is straightforward except for the presence of a much 
fainter star separated by only 1.5" on the sky. This was treated as 
third lig ht (Paper III). The mag nitude difference measured in the i 
band by IPaemgen et alj ( l200g) , M = 4.560 ± 0.017, was con- 
verted into a flux rat io to obtain the constraint on L3 . 

The data from IChan et alj (2011) cover five transits, two of 
them fully, with 765 datapoints. I fitted for Porb, and found that 
correlated noise was unimportant and that the LD-fit/fix solutions 
were good (Tab le A 128). These results were added to the B- and 
2-band ( Mandus hev et alj|2007h results from Paper III, and a com- 
paratively poor agreement was found. The most concordant param- 



© 0000 RAS, MNRAS 000, 000-000 



Studies of transiting extrasolar planets. V. 21 



1.00 



0.98 



E 0.96 



0.94 




. z 




-0.04 



-0.02 0.00 0.02 

Orbital phase 



0.04 



Figure 27. The B and z ^Ma ndushev et al.l |2007|) and i light curves 
iChan et al.l201lh of TrES-4 compared to the JKTEBOP best fits. Other com- 
ments are the same as Fig.|6] 



eter was, unusually, k {Xv = 0.4), whereas the other photomet- 
ric parameters have Xu ~ 1-1-1-8. The new i-band data yield re- 
sults roughly in the middle of those from the B and z data, which 
are formally discrepant. The final photometric parameters are the 
weighted average of the individual ones from the B, i and z data, 
with the errorbars inflated to enforce = 1 (Table A129). This 
leads to a significant difference with respect to Paper III, where the 
B data were rejected and the z-band results taken as final, but is in 
reasonable agreement with other studies. The best fits to the three 
light curves (i.e. those studied in Paper III as well as the current 
work) are plotted in Fig.l27l 

When calculating the physical properties of TrES-4 I found 
that the different model sets did not agree completely. The 
models are a strong outlier, but are also the closest to the 
spectroscopically-measured stellar T^s. These results have never- 
theless been rejected, and the final values based on those from the 
other four sets of models (Table A 130). They point to the planet 
having an unexpectedly large radius (i?b = 1.735 ± 0.072 ± 
0.007 Rjup), confirming the conclusions from previous studies. I 
recommend that new transit photometry and echelle spectra are 
obtained for TrES-4 in order to provide a definitive answer to the 
question of its physical properties. 

5.27 WASP-12 

At the time of its discovery, WASP-12 contained the hottest, 
most irradiated and shortest-period planet then known jHebb et al] 
l2009h . making it a key object in our understanding of gas-giant 
planets. Occultations of the planet by the star have been observed in 
the z band tLope z-Morales et al. 20_Ui), in the J, H and Ks bands 
dCroll et al.l2'oiih. and in t h e 3.6, 4.5, 5.8 and 8.0 /xm Sp itzer/IRAC 
bands jCampo et al.ll2oTTh . lMadhusudhan et all | |201 ih interpreted 
the occultation depths in these bands as evidence of an enhanced 
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Figure 28. Light curve of WASP-12 compared to the JKTEBOP bes t 
fits. From top to bottom the datasets are from iHebb et al] j2009l) . 
iMacieiewski et al.l t2011) . and two from IChan et alj j201lh . Other com- 
ments are the same as Fig.|6] 



ratio of carbon to oxygen in the atmosphere of the planet. Most 
recently, ICowan et alj ( 2012[) used Spitzer to obtain photometry at 



3.6 an d 4.5 /xm covering full orbits of WASP-1, and ISwain et al.l 
l l2012h have measured transit and occultation depths in the infrared 
using HSTAVFC 3. High-precision optical tran sit photome t ry has 
been obtained bv IMacieiewski et alj ( 1201 iT) and lChan et alj ( 1201 ih 
and used to confirm the large planetary radius measurement ob- 
tained bv lHebb et ai] (j2009). 

The possibility of orbital eccentricit y in th e WA SP-12 sys- 
tem has received much attention. .Hebb et al. I ( I2OO9.) found a 
value of e = 0.049 ± 0.015, which was unexpected because 
tidal effects should have circularised such a short-period orbit 
dGoldreich & Soteilll966l) . The prob ability that this ec centricity is 
real is 99.5% according to th e test of Lucv & Swe enev (1971). The 
occultation timing found by iLopez-Morales et"al. (2010.) leads to 
the constraint ecosctj = 0.0156 ± 0.0035 , which agrees with the 
eccentricity detection of'Hebb et al.' ('2009'). However, the later oc- 
cultation observations (Croll et al. 2011; Campo et al. 2011) define 
the eclipse shape much more precisely and constrain e to be very 
small and probably zero. Furth er support for a circu lar orbit comes 
from continued RV monitoring jHusnooetalj201lh . and from con- 
sideration of the effect of tidal deformation of the host star on its 
measured RV t Arras e t al. 2012). I have therefore assumed a circu- 
lar o rbit in the analysis below. 

iHebb et al] ( |2009|) presented photometry of two transits, a z- 
band one from LT/RATCam and a B-band one from the 0.8 m 
Tenagra telescope in Arizona. The former of the two datasets is 
of m uch better quality than the latter, and was included in my anal- 
ysis. |Maciei^s^eral](|20ll|) obtained photometry of two transits 
in the Johnson R band with the CAHA 2.2m/CAFOS. The second 
dataset suffers froin greater scatter, gaps and very little data outside 
transit, so I cons idered only the observations of the first transit. 
IChan et al] ( 120111) suffered from a lot of bad weather, but were able 



© 0000 RAS, MNRAS 000, 000-000 



22 John Southworth 



to get good observations of one transit in i using KeplerCarrF"! and 
one transit in V using NOT/ALFOSC. Botlr datasets were added to 
my analysis. 

One problem is the asphericity of the planet, which is techni- 
cally outside the range of validity of the EBOP model (models as- 
suming spherical planets will be more strongly affected). This was 
discussed in the context of OGLE-TR-56 in Sect. l5.2ll the corol- 
lary for WASP- 12 is that the volume-equivalent radius is up to 5% 
larger than the measured transit radius. 

For all datasets I found that the MC errors were bigger than 
the RP errors and that the LD-fit/fix solutions were good. Results 
for the 2-band data (613 points) are given in Table A131; for the R- 
band data (167 points) in Table A132; and for the i and V band data 
(470 and 670, r espectively) in Table A 133 and Table A 134. The R- 
band data from lMacieiewski et alj ( l201ll) are of very high quality, 
with a scatter of only 0.59 mmag per point, and for these it was 
possible to adopt the LD-fitted solution. The best fits are plotted in 
Fig.Hi 

Table A135 shows that the agreement between the four light 
curves is poor for k {Xv = 2.5) but fine for the other photometric 
parameters. Literature val ues are in good agr eement with my final 
results, except for those of lChan et alj j201 ih . These authors find a 
lower k than I achieve for any of the four datasets considered. 

The physical properties from the five different model sets dis- 
play a comparatively poor agreement. The VRSS models are the 
most discrepant, but not by enough to entertain rejecting them. My 
final results (Table A136) therefore have a significant systematic 
error component, which however is still dwarfed by the statisti- 
cal uncertainties. WASP-12 would really benefit from more pre- 
cise measurements of the stellar T^a and ; the mass of the star 
is uncertain by 14% (Ma = 1.38 ± 0.19 Mg) primarily due to 
the uncertainty in the star's atmospheric parameters. In compari- 
son to literature results, I find R^ = 1.825 ± 0.091 ± 0.024 R jup, 
which is 1(7 smaller th an that fromlMacieiewski et al. ll l20Tlh and 
la larger than that from lCham et al.l ( l201lh . WASP-12 b is therefore 
confirmed to have an unusually large radius. 
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Figure 29. The JGT (top) and LT (lower) light curves of WASP- 13 com- 
pared to the JKTEBOP best fits. Other comments are the same as Fig.|6] 



found photometric parameters which agree with my own, but with 
errorbars a factor of 2 smaller. 

The DSEP models returned results in some disagreement with 
the other model sets, so were discluded when calculating the final 
physical properties. The extant T^ff and [|^] measurements for 
WASP- 13 A are comparatively poor, and contribute to the uncer- 
tainties in the properties of this system. I find a larger and more 
massive star and planet, the change in each case being approxi- 
mately the same size as the errorbars (Table A140). My revised sys- 
tem parameters makes WASP-13b a rather large (i?b = 1.528 ± 
0.084 ± 0.004 Rjup) and low-density (p b = 0.134 ± 023 p jup) 
TEP. The co r respon ding values found bv lSkillen et al.l ( l2009h and 
iBarros et al. I l l2012l) are for a smaller and less massive star and 
planet. WASP- 13 deserves further photometry and spectroscopy. 



5.28 WASP-13 

ISkillen et alj ( |2009|) presented the discovery of the planetary na- 
ture of this system, with an analysis including observation of one 
full transit from the 0.95 m Jam es Gregory Telescope at St. An- 
drews, UK. iBarros et alj ( 1201 2|) subsequently observed four tran- 
sits using LT/RISE, although none in their entirety. I have analysed 
both datasets here, with Porb incuded as a fitted parameter for the 
LT data (1530 datapoints) but not for the single transit from the 
JGT (173 datapoints). In both cases I found that correlated noise 
was significant: the RP errors are greater than the MC errors by a 
factor of 1.5 for the JGT data and by a factor of 2 for the LT data. 
The full solutions can be found in Table AI37 and Table AOS. 

The photometric solutions are summarised in Table A 139 and 
plotted in Fig.[29] The agreement between the two datasets is good 
except for k, as usual. The parameters found bv lSkillen et al.l ^20091) 
disagree with my own, by roughly 2a depending on parameter. The 
reasons for this are not clear, given than I found the JGT and LT 
data to be consistent with each other, but may be related to inclu- 
sion of the SuperWASP data in their analysis. IBarros et al.l j20I2l) 

^" The observing log in lChan et al.l )201 ih specifies a z-band filter for this 
dataset but the datafile from ApJ labels it as i. J. Winn (2011, private com- 
munication) has confirmed that the latter is correct. 



5.29 WASP-14 

Discovered by Ijoshi et alj ( |2009|) . this planet is massive (7.9 ± 
0.5 Mjup) and has a definitely eccentric orbit (e = 0.087 ± 
0.002). These two prope rt ies are known to preferentially coincide 
jSouthworth et alj|20093) . I Johnson et alj ( [2009.) has found that the 
planetary orbit is misaligned with the stellar spin (A = —33.1 ± 
7.4 de grees) from Rossiter-McLaughlin observations. iBlecic et alj 
( I2OI l l) obtained Spitzer observations of three occutations in order 
to probe the atmospheric properties of the planet from its infrared 
flux distribution. 

The orbital eccentricit y was detected by [ joshi et alj ( |2009|) 
and has been confirmed by iHusnoo et alj ( I2OI1 I) from continued 
RV m onitoring. The most precise value comes from Blecic et al] 
( |20T Tj), whose three occultations allow a good constraint on 
ecoso;. I have adopted the measured ecoso; = —0.02557 ± 
0.00038 and e sin a; = 0.0831 ± 0.0021 as constraints on the JK- 
TEB OP models. 

IJosM et alj ( |2009l) presented an LT/RISE light curve of one 
transit containing 419 datapoints. The LD-fit/fix solutions are rea- 
sonable, and LD-fitted ones were not attempted (Table A141). Cor- 
related noise occurs, causing RP errors to be 1.5 times greater than 
the MC errors. 

iJohnson et aLl( l2009l) observed one transit (247 datapoints) us- 
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Figure 30. Light curves of WASP-14 from Ijoshi et alj j20O9l) and 
Ijohnson et al.l j2009l) compared to the JKTEBOP best fits. Other cormnents 
are the same as in Fig.|6] 



ing the UH 2.2m and OPTIC CCD, in a non-standard filter cen- 
tred on 850 nm and 40 nm wide. They used this to constrain the 
orbital ephemeris to aid their Rossiter-McLaughlin measurement, 
but did not attempt a revision of the other parameters of the sys- 
tem. I adopted LD coefficients for the Cousins / band as this is 
reasonably close to the actual filter used (given that the predicted 
LD coefficients are significantly scattered). I found that the MC er- 
rors were larger than the RP ones, and that the LD-fit/fix solution is 
good (Table A142). 

The light curve solutions are collected in Table A143 and the 
best fits are shown in Fig.|30] The agreement between the two light 
curves is unimpressive, with k having = 0.8 and the other 
photometric parameters having Xv ~ 1-5-1.9. There is no spe- 
cific reason to prefer one light curve over the other, so two results 
for the two have been combined and the errorbars inflated to ac- 
count for the modest discrepancy. WA SP-14 wo uld clearly benefit 
from further photometric observations. I Joshi et al. (2009) obtained 
results in agreement with those I find from their data, but in less 
good agreement with my final photometric parameters; no results 
are available from the other works cited above. 

It is no surprise that I find a significantly different set of 
physical properties of WASP-14 compared to previous works (Ta- 
ble A144). The star becomes more massive and more evolved 
(Ma = 1.35 ± 0.12 Mq and i?A = 1.67± 0.10 Rq). The planet's 
size is similarly revised, to Rh = 1.633 ± 0.092 ± 0.009 Rjup), 
an increase of 4o- over the value of 1.28llo;o82 Rjup found by 



|joshietalJ ( l200^ . This new i?b is in cle ar disagreement with stan- 
dard theoretical models of planets (e.g. iBodenheimer et aO l2003l : 
iFortnev et alj|2007l : lchabrier et al.ll2009h . New light curves are vi- 
tal in sorting out the moderate discrepancy between the two existing 
ones, and new stellar atmospheric parameters would allow a great 
improvement in our understanding of WASP-14. 



5.30 WASP-21 

iBouchv et all ( 1201 Ol) discovered this planetary system, and ob- 
tained its physical properties using their version of the dEB con- 
straint. They obtained adaptive-optics imaging which show no 
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Figure 31. The LT light curve of WASP-21 compared to the JKTEBOP best 
fit. Other comments are the same as Fig.|6] 



faint stars contaminating the flux from the system. The inter- 
est in WASP-21 lies in its Saturn-mass planet (Mb = 0.295 ± 
0.030 Mjup) and the metallicity which is the lowest known for a 
TEP host star ([§] = -0.46 ± 0.11). iBarros et al.l ( 1201 ih anal- 
ysed LT/RISE observations covering three transits, of which the 
first w as originally presented in the discovery paper i Bouchv et al.l 
I2OIOI) . They used the models to provide their additional con- 
str aint, yielding signifi cant differences compared to those found 
bv IBouchv etatl Ma = 0.86 ± O.O4M0 versus 1.01 ± 

0.03MQ,andMb = 0.27 ± 0.01 Mq versus 0.300 ± 0.011 Mq. 
WASP-21 is therefore an obvious candidate for inclusion in the Ho- 
mogeneous Studies project. 

The LT data contain 1302 points over three transits (one each 
in 2008, 2009 and 2010), of which only one has full coverage. 
These were modelled with Porb included as a fitted parameter. As 
with WASP- 13, 1 found the RP errors to be 1 .5 times larger than the 
MC errors, and that the LD-fit/fix solution was the best. The full set 
of solutions is in Table A145 and the best fit is plotted in Fig.|31| 
I find photometric par ameters in reasonable agreement with those 
of lBarros et al.l ( 12011 ) but with larger errorbars (Table A146). The 
parameters from Bouchv et al. ( 2010) are discrepant, as has been 
discussed by Barros et al. (201 1). 

The physical properties arising from the five model sets are 
not in great agreement and this is straightforwardly attributable to 
the low [§] , for which the model predictio ns diverge. I fin d stel- 
lar a nd planetary masses closer to those from lBouchv et al] ( I2OI0I) 
than lBarros et al.l ( 1201 ll) . but stellar and planetary radii which are 
bigger than found by either study (Table A 147). The planet is now 
one of the least dense known (pb = 0.137 ± 0.021 ± 0.003 p,7up). 
A new spectroscopic study of WASP-21 should be given priority, 
to confirm the low [^] and increase the precision of the Teff and 
[^] measurements. A new light curve would be useful, as the cur- 
rent one dominates the uncertainty in Ra and i?b • 



5.31 TEPs without new photometric studies: CoRoT-2, 
WASP-2, WASP-2, WASP-4, WASP-5, WASP-7, 
WASP-18 and XO-2 

CoRoT-2 was studied in Paper IV but has since been the subjec t 
of a high-resolution spectroscopic study bv lSchroter et alj ( I2OI lb . 
They find improved atmospheric parameters for the star {T^s = 
5598 ± 34 K and [^] = 0.04 ± 0.02). For comparison, the values 
used in Paper IV were 5696 ± 70 K and 0.03 ± 0.06, respectively. 
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The errorbars on the new measurements are rounded up to ±50 K 
and ±0.05 dex for the reasons given in Paper II. The new Tcff yields 
a slightly smaller system scale (Table AI48), which means lower 
masses and radii for both components. The Claret models disagree 
slightly with the others, resulting in comparatively large systematic 

errorbars for some calculated quantities. 

W ASP-1 is in receipt of new results from lAlbrecht et al.l 
(l201lh . who revised Ka to 125 ± 5ms"^ versus 111 ± 9ms"\ 
Teff to 6213 ± 51 K versus 6110 ± 50 K, and [§] to 0.17 ± 0.05 
versus 0.23 ± 0.08. It is slightly unusual that an upward revision 
in Tefi brings about a downward revision in [^] , as the two are 
positively correlated in the spectral analysis of F, G and K stars. 
The change in T^a is also more than the errorbars, so a confirma- 
tion of this result would be appropriate. Table AI49 shows that the 
new measurements modify the physical properties of WASP-1 only 
within the errorbars compared to Paper III. WASP-1 b's status as a 

hot and bloated planet is secure^ 

WASP-2 was studied bv lSouthworth et alj i2010h . New spec- 
troscopic results (cover ing both RVs and atm ospheric parameters) 
are now available from lAlbrecht et alj j20Ilh . It is not straightfor- 
ward to deduce the st ar's Teg as multiple values a vailable in the 
literatu re do not agre e : ICollier Cameron et al.l j2007l) f ound 5200 ± 
200 K ; iGhezzi et al.l feoioh obtained 5227 ± 4 6 K; iTriaud et all 
fcOlOh fo und 5150 ± 80 K;iAlbrecht et alj ( 1201 ih obtained 5206 ± 
50 K; and lMaxted eTal] j201lh deduced 5110 ± 60 K. I have there- 
fore adopted 5170 ± 60 K as a value which encompasses all these 
determinations and leans towards the slightly lower T^s from 
iMaxted et al.l j20Ilh . Their value wa s obtained usin g the InfraRed 
Flux Method (IRFM: Isiackwell & Shallis 1977; B lackwell et al] 
11980 *). which I consider to be the most reliable method due to its 
emphasis on the stellar spectral energy distribution (which is con- 
ceptually nicer for TeS determinations) and its limited dependence 
on theoretical calculations. TableA150 shows that the agreement 
between models is not great, as the star is substantially less massive 
than the Sun, but the overall results are a big step forward from pre- 
vious measurements. This is due to the more precise T^g and [^] , 
and means that WASP-2 is now one of the best-understood TEP 

systems. 

WASP-4 gets a new Tcfi of 5540 ± 55 K from lMaxted et al.l 



in the system parameters with the new Tcs. WASP-7 could still do 
with a better light curve and [|^] measurement. 

WASP-18 gets a new T^a of 6455 ± 70 K from lMaxted et al 



201 If), which is a substantial improve ment on the 5500 ± 100 K 



Gillon et alj2009l) used in Paper III and lSouthworth et alj j2009bh . 



The agreement between models is imperfect with respect to the sys- 
tem properties, but the final results achieve a high precision and are 
in agreement with literature values (Table AI5 1). WASP-4 is one of 
the best-understood TEP systems. The most effective improvement 
would be an increase in the precision of [^] , but this is not a high 
priority. 

WASP-5 has r e ceived a new Tcft measurement of 5770 ±65 K 
from lMaxted et all (l201ll) which is a useful improvement on the 
5700 ± 100 K i Gillon etal] |2009|) adopted in Paper III and by 
ISouthworth et al. I l2009al). The ensuing physical properties (Ta- 
ble AI52) show a modest agreement between models but a good 
agreement overall with literature values. WASP-5 could do with a 
better light curve and [^^] determination. New photometric data 
of WASP-5 have been obtained by several researchers - some is 
already published - and this will help in the near future. 

W ASP-7 has a revised Toft of 625 ± 70 K fromlMaxted et alj 
(|20II|). compared to the 6400± 100 K faellier et alj2009bl) usedTn 
ISouthworth et ai]l l20Ill) . Note that WASP-7 was recently found to 
have an almost pol ar orbit through obse rvations of the Rossiter- 
McLaughlin effect jAlbrecht etal]|20I2h . but this study did not 
generate a new if a. Table A153 shows that there is little change 



(201 which is more precise than the 6400 ± 100 K ( Hell ieretal 
i2009a) adopted in Paper III. Table A 1 54 shows that there is not 
much difference in the physical properties, but reinforces the sug- 
gestion that the system is young as a zero age is the preferred so- 
lution for three of the five models. The errorbars may be artifically 
small due to the concomitant edge effects. WASP-18 would bene- 
fit from a new spectral-synthesis study as well as additional light 

curve data. 

XO-2 has been awarded a revised Ka bv lNaritaetZI ( l201lh . 
who found 92.2 ± 1.7 ms~^ compared to the 85 ± 8ms~^ 
dBurke et al.l2007h used in Paper III. XO-2 A remains the TEP host 
star with the highest at ±0.45 ± 0.05 dex. The models con- 
tinue to disagree significantly (Table A155), due to the high [^] , 
and in such a way that is not accounted for by most analyses of this 
system. The new results are very similar to the previous ones from 
Paper III. 



6 PHYSICAL PROPERTIES OF THE TRANSITING 
EXTRASOLAR PLANETARY SYSTEMS 

The primary result of this work is measurements of the physical 
properties of 38 transiting planetary systems. These are collected in 
Table|3]for the stellar components and in Table|4|for the planetary 
components. When combined with previous results from the Ho- 
mogeneous Studies project, homogeneous physical properties are 
available for a total of 82 TEP systems. In the discussion below 
I have added results for all other known TEP systems as aggre- 
gated in the Transiting Extrasolar Planets Catalogue (TEPCa^j) 
on 2012/05/02, which comprises 199 TEPs and six transiting brown 
dwarf systems including those objects in the Homogeneous Studies 
project. 

Fig.|32|is a mass-radius plot for the host stars. This sample 
is very inhomogeneous due to observational biases and selection 
effects, particularly because the majority of objects have been dis- 
covered by wide-field surveys with small ground-based telescopes. 
The happy hunting ground is clearly unevolved stars in the mass 
interval 0.9-1.2 Mq. Lower-mass stars are intrinsically dim and 
predominate at faint apparent magnitudes, which are inaccessible 
to wide-field telescopes with coarse pixel scales. Higher-mass stars 
are brighter and plentiful in such surveys, but exhibit faster rota- 
tion and fewer spectral lines, making confirmatory RV measure- 
ments difficult. More evolved stars are easier for RV measurements, 
but their larger radii means transits of planet-sized objects are both 
longer and shallower, making it difficult to identify the photomet- 
ric signature of a TEP in the first place. Space-based surveys are 
less biased towards unevolved stars, so the characteristics of this 
diagram will change in the near future as new discoveries flood in 
from Kepler. 

Fig.|33] shows the mass-radius plot for the known TEPs in 
the region of parameter space occupied by Jupiter, with loci of 
constant density overlaid. Note that Jupiter does not lie on the 
Pb = /9.iup locus because of the rotationally-induced difference 
between its equatorial radius (7 1 492 km) and volume-equivalent 



The Transiting Extrasolar Planet Catalogue (TEPCat) is available at 
ihttp : / /www .astro ■ keele ■ ac ■ uk/ jkt/tepcat/ 
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Table 3. Derived physical properties of the stellar components of the TEPs studied in this work. For most quantities the first eiTorbar gives the statistical errors 
and the second errorbar gives the systematic errors arising from the dependence on stellar theory. 



System 



Semimajor axis (AU) 



Mass (Mq) 



Radius (R0) 



log gA [cm/s] 



Density (pq) Age (Gyr) 



CoRoT-2 

CoRoT-17 

CoRoT-18 

CoRoT-19 

CoRoT-20 

CoRoT-23 

HAT-P-3 

HAT-P-6 

HAT-P-9 

HAT-P-14 

Kepler-7 

Kepler- 12 

Kepler-14 

Kepler-15 

Kepler-17 

KOI- 135 

KOl-196 

KOl-204 

KOl-254 

KOI-423 

KOI-428 

OGLE-TR-56 

OGLE-TR-111 

OGLE-TR-113 

OGLE-TR-132 

OGLE-TR-L9 

TrES-4 



0.02835 ± 0.00029 
n n/i8in +"-00109 

0.U4S1U _o.oo064 

0.02860 ± 0.00065 
u.uoi^ _o.ooii 
0.0892 ±0.0028 

0-04809 

0.03842 ± 0.00050 
0.05244 ± 0.00077 
0.0529 ±0.0014 
0.06108 ±0.00069 
0.06324 ± 0.00042 

U.UOOO _o 0013 

0.0771 ±0.0010 
0.0583 ±0.0022 
0.02605 



+0.00018 
0.00124 

0.0440 ±0.0013 

n n9SQc; +0.00091 

U.U2Sy5 _o.00086 

0.0457 ±0.0017 
0.02954 ± 0.00068 
0.1539 ±0.0073 
0.0802 ±0.0032 
0.02453 ± 0.00046 
0.04676 ± 0.00092 
0.02265 ±0.00099 
0.03026 ± 0.00098 
0.04047 ± 0.00093 
0.0502 ±0.0010 



WASP- 


-1 


osq20 +0-00053 
u.uaazu „o. 00031 


WASP- 


-2 


0.03092 ± 0.00046 


WASP- 


-4 


0.02318 ±0.00037 


WASP- 


-5 


0.02740 ± 0.00040 


WASP- 


-7 


0.0624 ±0.0011 


WASP- 


-12 


0.02309 ± 0.00096 


WASP- 


-13 


0.0557 ±0.0018 


WASP- 


-14 


0.0372 ±0.0011 


WASP- 


-18 


0.02043 ± 0.00029 


WASP- 


-21 


0.0516 ±0.0020 



XO-2 



03616+°°°°'''^ 
U.UODiO _o. 00182 



± 0.00037 

+0.00023 
-0.00034 

±0.00019 

+0.0003 
-0.0003 

±0.0004 

+0.00039 
-0.00056 

± 0.00063 
±0.00051 
± 0.0003 
± 0.00036 
± 0.00090 

+0.0003 
-0.0002 

± 0.0006 
±0.0009 

+0.00037 
-0.00050 

±0.0002 

+0.00020 

-0.00016 

±0.0001 
± 0.00075 
±0.0013 
± 0.0032 
±0.00041 
±0.00081 
±0.00021 
± 0.00029 
± 0.00026 
±0.0002 

+0.00015 
-0.00015 

± 0.00038 
±0.00028 
±0.00019 
± 0.0003 
±0.00030 
±0.0001 
± 0.0002 
±0.00014 

±0.0012 

+0.00100 

-0.00167 



0.997 ±0.030 

1 043+0.073 
1.U10 _o.041 

0.861 ±0.059 

, ,01 +0.125 
i-J-Ol _o.074 

1.11 ±0.10 

-, -190+0.088 
1-1^^ -0.061 

0.900 ±0.036 
1.295 ±0.057 

1.28 ±0.10 
1.418 ±0.048 

1.413±0.081 

1 -|.;7 +0.147 
1.101 -0.078 

1.318±0.052 
1.08 ±0.12 

l.uou _o.i5o 
1.24 ±0.11 

n Qqq +0-092 

u.yjy _o.o82 

1.21 ±0.13 
0.570 ± 0.040 
1.08 ±0.16 
1.46 ±0.18 
1.339 ±0.077 
0.846 ± 0.049 
0.755 ±0.095 

1.29 ±0.12 
1.427 ±0.099 
1.339 ±0.084 

1 +0.052 
l.ZDO „o.030 

0.851 ±0.038 
0.927 ± 0.044 
1.033 ±0.045 
1.317 ±0.069 
1.38 ±0.18 

1.22 ±0.12 
1.35 ±0.12 
1.274 ±0.054 
0.98 ±0.12 

024+0-054 
u.a^i -0.122 



± 0.039 

0.015 
0.022 
±0.017 
+0.024 
-0.022 
±0.01 
+0.027 
-0.039 

± 0.044 
±0.038 
±0.02 
± 0.025 

± 0.060 

+0.022 
-0.011 
± 0.029 
±0.05 

+0.045 
-0.061 

±0.02 

+0.019 
-0.015 

±0.01 
± 0.044 
±0.03 
±0.18 
± 0.066 
± 0.043 
±0.021 
±0.04 
± 0.028 
±0.016 

-0.015 
-0.014 
± 0.032 

± 0.034 

±0.021 

± 0.022 

±0.05 

±0.01 

±0.02 

± 0.027 

±0.07 

+0.076 
-0.124 



0.901 ±0.015 ±0.012 

-1 cn +0.47 +0.01 
-^•"^ -0.13 -0.01 

0.924 ± 0.057 ± 0.006 

1 c;7C +0.354 +0.011 
l.OlU _o.o96 -0.010 

1.34 ±0.37 ±0.01 

1 +0.26 +0.01 
-^•'^ -0.16 -0.02 

0.870 ±0.016 ±0.014 
1.518±0.070±0.015 
1.339 ±0.080 ±0.008 
1.591 ±0.056 ±0.009 
2.028 ± 0.025 ± 0.029 
1 AQn +0.061 +0.009 
l.^iau -0.038 -0.005 

2.09 ±0.11 ±0.02 
1.253 ±0.047 ±0.020 

n QSq +0-013 +0.014 
U.MOO _o.o47 -0.019 

1.332 ±0.046 ±0.007 

n Q^Q +0.062 +0.006 
U.yoa _o.031 -0.005 

1.46 ±0.12 ±0.00 
0.539 ±0.036 ±0.014 
1.225 ±0.092 ±0.010 
2.48 ±0.17 ±0.10 
1.737 ±0.035 ±0.029 
0.822 ±0.019 ±0.014 

0. 766 ±0.038 ±0.007 
1.338 ±0.075 ±0.013 
1.499 ±0.042 ±0.010 
1.834 ±0.087 ±0.007 

1 AR<^ +0.051 +0.006 

1. 'iDO -0.079 -0.005 

0.823 ±0.015 ±0.010 
0.910 ±0.015 ±0.011 
1.088 ±0.040 ±0.008 
1.478 ±0.088 ±0.008 
1.619 ± 0.076 ±0.021 
1.657 ± 0.079 ±0.004 
1.666 ±0.097 ±0.009 
1.228 ±0.042 ±0.009 
1.186 ±0.081 ±0.028 

n QRO +0.026 +0.027 
U.yoz _o.054 -0.044 



4.527 ±0.1 

4.035 to 
4.442 ± 0.1 
4.115 1« 
4.23 ±0. 

4.008^0 
4.513 ±0.1 
4.188 ±0.1 

4.293 ± 0. 
4.187 ±0. 
3.974 ± 0.1 
4.155 
3.918 ±0.1 
4.275 ± 0. 
4.481 
4.281 ±0. 
4.466 
4.192 ±0.1 
4.730 ±0.' 

4.294 ± 0. 
3.812 ±0. 
4.086 ±0.1 

4.536 ±0.1 
4.548 ± 0.1 
4.297 ± 0.1 
4.241 ±0. 
4.038 ±0.1 
4.209 

4.537 ±0.1 
4.487 ±0.1 
4.379 ± 0. 
4.218 ±0. 
4.159 ±0. 
4.086 ±0.1 
4.126 ±0.1 
4.365 ± 0. 
4.281 ±0. 
4.436 



1+0 



016 ±0.006 

.069 +0.002 
208 -0.003 

043 ± 0.003 

.047 +0.003 
.152 -0.003 

24 ±0.00 
.076 +0.004 
.107 -0.005 

020 ± 0.007 

035 ± 0.004 
046 ± 0.002 

025 ± 0.003 
018 ±0.006 

.018 +0.003 
.014 -0.001 

040 ±0.003 

017 ±0.007 

.005 +0.006 
.023 -0.008 

017 ±0.002 
.018 +0.003 
.054 -0.002 

055 ±0.001 
039 ±0.011 
039 ± 0.004 
048 ±0.017 

018 ±0.007 
024 ± 0.008 

021 ±0.004 

036 ±0.004 
016 ±0.003 
033 ± 0.002 

.051 +0.002 
028 -0.002 

016 ±0.005 
009 ± 0.005 

030 ± 0.003 
048 ± 0.002 
023 ± 0.006 
033 ±0.001 
042 ± 0.002 
027 ± 0.003 

031 ±0.010 

.037 +0.012 



0.243^ 



1.362 ±0.064 

-0.066 
-0.126 

1.09 ±0.16 

U.OUZ_o 125 

0.46 ± 0.48 

214+° °®'' 
^■^l'*_0.066 

1.365 ±0.078 

0.370 ± 0.045 

0.534 ± 0.082 

0.352 ±0.031 

0.169 ±0.007 

3496+° °°®'' 
u.o'iau_o 0091 

0.145 ±0.020 

0.548 ± 0.008 

-I -19-1+0.015 
1-1^1-0.034 

0.523 ±0.022 

-I ,„^+0.046 
1-1O0-0.183 

0.388 ±0.071 
3.63 ±0.52 
0.586 ± 0.072 
0.096 ±0.014 
0.256 ±0.011 
1.52 ±0.10 
1.679 ±0.064 
0.540 ± 0.062 
0.424 ±0.018 
0.217 ±0.025 

n 4n-^+o-o69 

U.4Uci_o 037 



1, 


,524 


± 


0, 


,067 


1, 


,230 


± 


0, 


,022 


0, 


,801 


± 


0, 


,080 


0, 


,408 


± 


0, 


,068 


0, 


,325 


± 


0, 


,016 


0, 


,268 


± 


0, 


,029 


0, 


,293 


± 


0, 


,042 


0, 


,687 


± 


0, 


,062 


0, 


,587 


± 


0, 


,061 



,+1.8 +0.7 
-1.7 -1.2 
,+0.7 +0.3 
-3.5 -0.4 
unconstrained 
. g+0.6 +0.3 
-1.7 -0.3 
fl.6 +0.8 
-11.4 -1.6 
+2.4 +0.5 
-1.7 -1.0 
+4.2 +3.6 
-3.8 -2.7 
+0.4 +0.4 
-0.8 -0.3 
-1 (.+1.4+0.4 
-2.0 -0.4 
-I +0.4 +0.4 
1-'-' -0.4 -0.2 
+0.0 +0.3 
0.1 -0.3 

1.5 +0.3 
5.0 -0.3 

1.6 +0.9 



1-2: 

8.0^ 



5.9; 

5.4 
7.5 
2.3 



2.9; 
5.3" 



6.2: 
4.7" 



+0 

-0.035 -0.020 



1 034+°-l2'' 
l.uo^_0 058 



-1.9 -0.5 
f, 7+6.7 +0.9 
' -1.6 -1.1 
-I c: +9.9 +2.6 
-^■^ -1.2 -1.5 
r, 7+2.4 +1.2 
^- ' -1.4 -0.5 
+5.2 +0.6 
■5.1 -0.5 
+4.4 +0.3 
■6.7 -0.4 
unconstrained 
4 4+2.2+0.7 
^■^ -4.2 -0.5 
n s+3-6 +0-7 
U-»_i.5 _o.8 
q n+2-9 +0.3 
•^■^ -1.3 -0.4 
Q 7+9.6 +5.3 
^- ' -8.8 -2.8 
unconstrained 

1.7: 

i.o: 
2.5: 
2.7: 
8.2: 
6.4: 
5.6: 
1.9: 
2.0: 
5.0: 
2.4: 
0.4: 

unconstrained 

(, c +23.2 +5.3 



,5 


+2, 


1 


,2 


-0, 


,7 


,6 


+0, 


,3 


,3 


-0, 


,2 


,8 


+0, 


,9 


,7 


-2 


,5 


,4 


+0. 


,5 


,5 


-0, 


,4 


,3 


+2, 


,5 


,3 


-1 


,7 


,7 


+ 1. 


,9 


,3 


-1 


,9 


,2 


+ 1, 


1 


,2 


-1 


,0 


,7 


+0, 


,3 


,6 


-0, 


,3 


,7 


+0, 


,3 


,6 


-0, 


,6 


,6 


+1. 


,6 


,7 


-0, 


,9 


,5 


+0, 


2 


,0 


-0, 


,2 


,0 


+0, 


4 


,7 


-0, 


,3 



-3.5 



-3.8 



radius (66 991 krrF^. This Figure shows that the dominant popu- 
lation of TEPs is currently of mass 0.5-1.2 Mjup and radius 0.9- 
l.VRjup. The Satumian objects are currently easily differentiated 
based on their masses and radii, but future discoveries may fill in 
the gaps between them, the Jupiter-like planets, and the lower-mass 
planets. 

Correlations have prev iously been notice d between Pqa an d 
gb jSouthworth et al.ll2007 ll and Porb and Afb jMazeh et ai]|2005b . 
The relevant plots are shown in Figs.[34]and|35] In both cases there 
are quite a few planets whose properties put them outside the range 
of -Porb shown in these plots. For the following discussion I have 
neglected those planets which lie outside the relevant plot due to 
large mass or long period. I have also ignored super-Earth planets 



^'^ |http ://nssdc.gsfc.nasa. gov/planetary /fact sheet/] 
jupiter f act . html 
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(Mb < 10 M^) in order to obtain a sample of definitely gaseous 
objects. 

The rank correlation test o f lSpearmanl ( ll904h returns a prob- 
ability of 99.988% (3.9a) that the Porb-ffb correlation is real and 
99.987% (3.8cr) that the Porb-^^b correlation is real. In both cases, 
inclusion of all TEPs returns an increased significance level of 
the correlation (4.8(t and 5.8a" respectively). Both correlations are 
therefore supported by the current data. 

[Hansen & Barmaril ( l2007h divided up eighteen of the twenty 
TEPs then known into two classes based on their position in a di- 
agram of O versus Tcq. An updated version of the diagram can 
be seen in Fig.[36l and agrees with previous conclusions (Paper II) 
that the division between the classes is blurred into insignificance. 
A dotted line at O = 0.055 has been drawn to show the ex- 
pected boundaries between Class I (O ~ 0.07 ± 0.01) and Class II 
{Q ~ 0.04 ± 0.01). The previously postulated gap in the distribu- 
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Table 4. Derived physical properties of the planets of the TEPs studied in this work. For many quantities the first eiTorbar gives the statistical eiTors and the 
second errorbar gives the systematic eiTors arising from the dependence on stellar theory. 



System 


Mass(Mjup) 


Radius ( R,jup) 


9b (ms 2) 


Density (p,jup) 


T^c'q (K) 






CoRoT-2 

CoRoT-17 

CoRoT-18 

CoRoT-19 

CoRoT-20 

CoRoT-23 

HAT-P-3 

HAT-P-6 

HAT-P-9 

HAT-P-14 

Kepler-7 

Kepler- 12 

Kepler- 14 

Kepler- 15 

Kepler-17 

KOI- 135 

KOl-196 

KOl-204 

KOl-254 

KOI-423 

KOI-428 

OGLE-TR-56 

OGLE-TR-111 

OGLE-TR-113 

OGLE-TR-132 

0GLE-TR-L9 

TrES-4 

WASP-1 

WASP-2 

WASP-4 

WASP-5 

WASP-7 

WASP- 12 

WASP- 13 

WASP-14 

WASP- 18 

WASP-2 1 

XO-2 



3.06' 



3.57 ±0.13 

^•^O -0.24 

3.27 ±0.17 
1 090 

I.UJU _Q 077 

5.06 ±0.36 

. +0.32 
-0.30 

0.584 ±0.020 
1.063 ±0.053 
0.778 ±0.083 
2.271 ±0.079 
0.453 ± 0.067 
430 +°"^^ 



7.68 ±0.37 
0.696 ± 0.097 
2.340 

3.09 ±0.21 

u.tjo _o.070 

1.030 ±0.087 
0.500 ±0.055 
17.9 ±1.8 
±0.40 
±0.17 
±0.10 
±0.20 
±0.15 
±1.5 
0.897 ±0.075 
0.980 

0.880 ±0.031 
1.249 ±0.043 
1.595 ±0.048 
0.98 ±0.13 
1.43 ±0.13 
0.512 ±0.060 
7.90 ±0.46 
10.38 ±0.30 
0.295 ± 0.027 

593+0 025 



2.16 
1.41 
0.55 
1.23 
1.17 
4.4 



±0.09 

+0.02 
-0.04 

±0.04 

+0.015 
-0.014 

±0.04 

+0.05 
-0.07 

±0.019 
±0.021 
± 0.009 
± 0.027 
±0.013 

+0.005 
-0.003 
±0.11 

±0.022 

+0.066 
-0.089 
±0.03 
+0.007 
-0.005 

±0.006 
±0.026 
±0.3 
±0.18 
±0.05 
±0.02 
±0.02 
±0.02 
±0.1 
± 0.007 

+0.008 
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tion is no longer apparent, so the division into t-wo classes can now 
be safely ignored. 

Fig.[37] takes a different approach to illustrating the distribu- 
tions of physical properties of the known TEPs. Distributions have 
been obtained for each property given in Table|4l and for the Homo- 
geneous Studies objects and for all objects. Each distribution was 
obtained by calculating a Gaussian function for each planet with a 
mean and standard deviation equal to the relevant property and its 
uncertainty, and then taking the sum over all planets in the sample. 
This produces an alternative to the histogram which avoids the loss 
of information implicit in the binning process. 

Fig.l37b shows that the mass distribution of TEPs has a broad 
plateau around 0.5-1.2 Mjup, with a peak at the lower edge of this 
interval. A secondary peak, representing the Saturnian TEPs, oc- 
curs around 0.2 Mjup. The distribution then rises towards the low- 
est masses. Fig.l37b shows the situation for planet radius; the dis- 
tribution is bimodal with peaks near 1.2Rjup and 0.2Rjup. The 
lower-radius peak is almost entirely absent from the Homogeneous 
Studies sample, as it is dominated by planets from the Kepler satel- 



lite which were discovered and characterised by TTVs. Such sys- 
tems are more complicated to study than normal Hot Jupiters, and 
there are no plans to include them in the current series of papers. 

Fig.l37b shows the distribution of planet surface gravities. 
The sharp peak around 9ms^2 (jug to the extremely precise 
measurement for HD 209458. The determination of g\, requires 
only directly measurable photometric and spectroscopic parame- 
ters (rb, i, Ka and e), so can be measured to extremely high 
precision in favourable cases. In the case of HD 209458, gb = 
9.30 ± 0.08 ms^2 (Paper I). The distribution of planetary densi- 
ties is plotted in Fig.l37tl. and has a very similar functional form to 
that for gb but with greater smoothness. This is because the pre- 
cision of a pb value is very sensitive to the quality of the tran- 
sit light curves, but is also subject to systematic error from the 
model-derived additional constraint. In the case of HD 209458, 
Pb = 0.254 ± 0.004 ± 0.002 p,jup. 

Fig.|37b depicts the distribution of T^^, which is noisy but has 
a broad peak around 1500-1800 K. Ground-based surveys are ob- 
servationally biased in favour of short-period planets, so hot planets 
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Figure 32. Mass-radius plot for the host stars of the known transiting ex- 
trasolar planets. Those objects studied in this work are shown with (blue) 
filled circles and numbers taken from the literature with (red) crosses. 




Figure 33. Mass-radius plot for the known transiting extrasolar planets. 
Symbols are as in Fig. |32| The four gas giant planets in our Solar System 
are denoted with (black) filled circles. Dotted lines show loci where density 
equals 1.0, 0.5 and 0.25 p.jup. 



are over-represented in the known TEP population. Finally, Fig.|37f 
shows the result for The division between Class I and Class II 
planets jHansen & Barmanll2007h is shown with a grey dotted line, 
and coincides with one of several local minima in the distribution 
function. 



7 FOLLOW-UP OBSERVATIONS 

Most of the TEPs in the current work would benefit from further 
observations of some sort, and this is summarised in Table[5] In 
many cases the dominant uncertainty stems from the quality of the 
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Figure 34. Plot of the orbital periods versus the surface gravities of the 
dominant population of known TEPs. Symbols are as in Fig. |32| 
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Figure 35. Plot of the orbital periods versus the masses of the dominant 
population of known TEPs. Symbols are as in Fig. 1321 



light curve. This remains true for many of the CoRoT systems, even 
though they have space-based data. In each case this is due to their 
relative faintness (remember that the aperture of the CoRoT tele- 
scope is only 27 cm) and/or that few transits were seen due to either 
a short observing sequence or a long orbital period. 

Additional RV measurements are useful too. In many cir- 
cumstances, particularly for the fainter objects, the RVs are good 
enough to unambiguously confirm the planetary nature of a system 
but are the dominant source of uncertainty in the planetary masses. 
Now over 200 TEPs are known it seems appropriate to concen- 
trate follow-up resources on measuring the physical properties of a 
golden subset of these to high precision. An additional requirement 
of RVs is definition of the orbital shape (e and uj), and imprecise 
measurements of these quantities compromise measurements of the 
photometric parameters (in particular ta). 
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Figure 36. Plot of equilibrium temperature versus Safronov number for the 
full sample of planets. Objects shown with (blue) circles were studied in this 
work; those which are just (red) errorbars were not. Th e dotted line shows 
the sep aration between Class A and Class B proposed bv lHansen & BarmanI 
\200i). 



The physical properties of quite a few of the TEPs are also 
limited by the precision of the Tea and measurements avail- 
able. In many cases this can be improved, but in some cases this is 
not an option because the errorbars are already close to the limit set 
by our understanding of low-mass stars (taken to be 50 K in Teg 
and 0.05 dex in [^]). There is no immediate prospect of lower- 
ing these threshold s; in fact there is evidence that they are already 
slightly optimistic ( Bruntt et al]|2010bll2012h . 

Quite a few of the TEPs have ephemerides which will become 
uncertain over the timescale of a few years. To investigate this I 
compiled a catalogue of ephemerides of all known TEPs and iden- 
tified the first predicted times of transit which were uncertain by 
one hour, and by half of one transit duration. A list of objects for 
which one of these dates is earlier than the year 2022 is given in Ta- 
ble|6] The Kepler planets are not included in this analysis because 
they continue to be observed and will benefit from significantly im- 
proved ephemerides once these data become available. The list of 
the other planets is dominated by CoRoT objects, and it is notable 
that the ephemerides for CoRoT-4, CoRoT-14 and CoRoT-17 are 
already uncertain by more than one hour. HAT-P-31 is present in 
the list due to the fact that there is no published follow-up pho- 
tometry of this object at all, only the discovery data from the HAT 
survey cameras. Further photometric observation s of these objects , 
such as those obtaine d by the TERMS project teane et alj|2009l : 
IPragomir et al.ll201ll) . are needed before the predictive power of 
their ephemerides deteriorates much further. 



8 DISCOVERY RATE AND SKY POSITIONS OF THE 
TRANSITING EXTRASOLAR PLANETARY SYSTEMS 

Whilst scientifically less relevant than the analyses described 
above, the discovery rate of the transiting planets is of sufficient 
interest to be worthy of discussion. Fig.[38] shows the discovery 
rate of the known TEPs. For discovery time I have taken the date 




Figure 37. Plots of the distributions of planet properties. From top to bottom 
these are M\,, R\^, g^,, p^. T^c'q ™d 0. In each case the upper (red) line 
shows the distribution for all known TEPs and the lower (blue) line shows 
that for the TEPs considered within the current series of papers. 



at which important detail^j of that system became available in a 
refereed journal article. If the article was lodged on the arXi\EI 
preprint server prior to this, I have instead taken the date at which 
this preprint became public. 

The population of known TEPs continues to increase at an 
exponential rate; a plot of Fig.[38]but with a logarithmic ordinate 
axis (not shown) exhibits an approximately straight line. The early 
front-runners were the OGLE survey (including the follow-up ob- 
servations which were performed by a range of assoc iated and inde- 
pende nt consortia), but the SuperWASP consortium jPoUacco et al.l 
I2OO6I) held the lead during the years 2009-1 1. The greatest number 
of known TEPs is now attributable to data obtained by the Kepler 
satellite, and once again this number includes both associated and 
independent groups of researc hers. Kepler has no w been used to 
find over 2000 candidate TEPs jSatalha et alj2012h and possible as 



^■^ I define the 'important details' to comprise confirmation of the planetary 
nature of the system (which in turn requires an assessment of the physical 
properties of the putative planet) and sufficient information to allow the 
execution of follow-up observations (i.e. the full sky position and orbital 
ephemerides). 

,http : / /arxiv . org/| 
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Figure 38. Plot of the discovery rate of the known TEPs. The total number is shown by the upper black hne, whereas the discovery rate of selected individual 
surveys are indicated below this using coloured lines as indicated in the key on the top left of the plot (please see the electronic version of this paper for the 
colour figure). 



many as 5000 dXenenbaum et al.l2012h : the rate of false positives is 
probably small but remain s under discussion jMorton & Johnsoril 
1201 ll : ISanteme et al.ll2012l : IColon et alj|2012l) . The candidate Ust 
provides the ammunition to continue the exponential rise in the 
number of known TEPs. 

Fig.|39l shows the sky positions of the known TEPs, overlaid 
with the celestial equator and galactic plane, and is a good tool to 
illustrate the strength and diversity of the selection effects which 
afflict this population. The most obvious feature in Fig.[39] is a 
strong clustering around 19hr right ascension and +42° declina- 
tion, which is the Kepler field. The Kepler satellite jBorucki et al.l 
I2OIOI) is currently the most successful TEP discovery machine, and 
has a modest field of view (105 deg^). This concentration of points 
is the result. 

Additional clusters of TEPs occur around the intersections of 
the Galactic plane with the ce lestial eq uator, the two areas tar- 
geted by the CoRoT satellite ( Baglin et al]|2006 ). CoRoT aims for 
these positions in order to obtain a high surface density of stars 
(to give plenty of targets in its small 8.2 deg^ field of view) which 
are easily followed up using ground-based telescopes in both hemi- 
sph eres. Two other smal l groups of TEPs come via the OGLE sur- 
vey jUdalski et al.l2002cl) observations of the Galactic centre ( 18 hr 
R.A. and —30° Dec) and Carina regions (llhr R.A. and —60° 
Dec). 

The wider spread of objects in Fig.[39]is dominated by HAT 

teakos et al.. ,2002) in the north and SuperWASP tPoUacco et all 



1 20061) in the soutlF^. These, plus other similar projects, operate 
multiple small wide-field telescopes with very coarse pixel scales. 
They are thus able to cover large amounts of sky but must avoid the 
crowded fields towards the Galactic plane, due to problems with 
blending. Blending makes it more difficult to measure good pho- 
tometry and also elevates the rate of false alarms in a search for 
TEPs. 

The surveys operating small ground-based telescopes tend 
to point away from the Galactic plane, which introduces a bias 
towards intrinsically faint stars. This is advantageous because it 
reduces the fraction of giants in the studied sample; these sur- 
veys cannot find planets around giant stars but are able to detect 
larger transiting object which are of lesser scientific interest. Con- 
versely, the Kepler satellite studies stars with fainter apparent mag- 
nitudes which are located close to the Galactic plane. These stars 
have slightly different population charact eristics to those studied 
by other surveys with different designs; Schlaufman & Laughlm| 
( 201 1) show that the host stars of the Kepler planet candidates are 
preferentially metal-rich. Such biases must be accounted for when 
studying the population characteristics of the known TEPs and their 
host stars. 



^ Declaration of bias: the author is a member of the SuperWASP consor- 
tium. 
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Figure 39. Plot of the sky positions of the known TEP systems. The celestial equator is shown using a solid grey line and the plane of the Milky Way galaxy 
with a dashed grey line. The TEPs are colour-coded to indicate those found by selected discovery surveys (please see the electronic version of this paper for 
the colour figure). 



9 SUMMARY 

The Homogeneous Studies project aims to provide measurements 
of the physical properties of a large sample of transiting extraso- 
lar planetary systems, using consistent methods and with careful 
attention paid to the estimation of robust statistical and systematic 
errors. The transit light curves are modelled using the JKTEBOP 
code, and uncertainties are gauged with Monte Carlo and residual- 
permutation algorithms. Attention is paid to the treatment of limb 
darkening, orbital eccentricity, and contaminating 'third' light. To 
the resulting photometric parameters are added measured spectro- 
scopic quantities: the Tea, and velocity amplitude of the host 
star. One additional constraint is needed, and is supplied either by 
one of five sets of theoretical stellar evolutionary models or by a 
semi-empirical calibration of low-mass star properties based on de- 
tached eclipsing binaries. The statistical errors are propagated us- 
ing a perturbation algorithm, leading to complete error budgets for 
each output quantity. The inclusion of multiple stellar models al- 
lows systematic errors to be deduced too. 

The current paper has presented complete analyses of thirty 
TEP systems, and updated the physical properties of eight more 
systems based on newly published spectroscopic results. Combined 
with previous work, this yields a total of 82 transiting planets and 
host stars with homogeneously-measured physical properties. In 
many cases these results are based on more numerous datasets than 
previous measurements, so should be preferred over literature val- 



ues even if homogeneity is not of specific importance to the matter 
in hand. Headline results from the current paper are summarised 
below. 

Analyses in the current work were the first to consider all 
available good light curves for HAT-P-3 (five datasets), HAT-P- 
6 (three), HAT-P-9 (four), HAT-P-14 (four), WASP- 12 (four) and 
WASP- 14 (four). For WASP- 14 this has resulted in a change in 
planetary radius (7?b) from 1.28lt|5 [;8^ Rjup to 1.633 =E 0.092 zE 
0.009 R,iup, moving this object away from the parameter space oc- 
cupied by theoretical predictions. 

The current study is also the first to present results based 
on short-cadence data for Kepler-14, Kepler-15 and KOI-135. The 
measured mass of Kepler-14 A has decreased by 2.6(t, and the mea- 
sured radii of both components in the Kepler-15 system are both 
significantly larger. The density deduced for Kepler-15 b is lower 
by a factor of two compared to previous studies. 

Additional data was modelled for KOI- 196, KOI-204, KOI- 
254, KOI-423 and KOI-428. In two cases (KOI-423 and KOI-428) 
the available data is more than a factor three more extensive than 
used in published studies. For KOI-428 the radii of the compo- 
nents are larger than previously found; the star is now comfortably 
the largest known known to host a TEP at Ra = 2.48 =E 0.17 =E 
0.20 Rq. 

The physical properties of OGLE-TR-56 have been heavily 
revised based on the first complete analysis of newly-published 
high-quality photometry cAdams et al. ,201 lb) . The radius of the 
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Table 5. Summary of which further observations would be useful for the 
TEPs studied in this work. * indicates where additional data would be use- 
ful, and denotes where it would be useful but difficult to either obtain or 
analyse. I do not suggest the need for new data when such data are already 
in the process of being obtained (e.g. Kepler photometry). 



System 


Photometric 


Kaclial 


Spectral 




observations 


velocities 


synthesis 


CoRoT-2 








CoRoT-17 




** 


** 


CoRoT-18 


** 






CoRoT-19 


* 






CoRoT-20 


* 






CoRoT-23 


** 


** 


** 


HAT-P-3 








HAT-P-6 


* 






HAT-P-9 


* 




* 


HAT-P-14 


* 






Kepler-7 




★ 




Kepler- 12 








Kepler- 14 








Kepler- 15 




■* 


■* 


Kepler- 17 








KOI- 135 






■*• 


KOI- 196 






■* 


KOI-204 




★* 


★* 


KOI-254 




★* 


★* 


KOI-423 


* 


■* 


★ 


KOI-428 




** 


** 


OGLE-TR-56 




★* 




OGLE-TR-111 




★* 




OGLE-TR-113 




★* 


★* 


OGLE-TR-132 




■** 




OGLE-TR-L9 




** 


** 


TrES-4 




■*• 


■*• 


WASP-1 






■*• 


WASP-2 








WASP-4 








WASP-5 






★ 


WASP-7 


* 


* 


■*• 


WASP- 12 






■*• 


WASP- 13 




★ 


■* 


WASP-14 






■*• 


WASP- 18 








WASP-2 1 


* 




★ 



Table 6. Limits of the current ephemerides of the known TEPs. The two 
dates for each TEP indicate the first transits whose midpoints are uncertain 
by 1 hour and by half the transit duration. RHJD = HID - 24000000. The 
Kepler planets are not included because they are the subject of continued 
observations using the satellite. 



TEP 



1 hour uncertainty 
RHJD UT date 



Half-transit uncertainty 
RHJD UT date 



CoRoT-04 

CoRoT-14 

CoRoT-17 

CoRoT-20 

HAT-P-3 1 

CoRoT-23 

CoRoT-10 

CoRoT-09 

OGLE-211 

CoRoT-08 

CoRoT-15 

Qatar- 1 



55181.195810 
55268.529920 
55447.075200 
56550.756250 
56573.327850 
56823.191100 
57040.629000 
57461.558700 
58524.988640 
58780.283621 
59004.400840 
59217.596165 



2009 12 15 
201003 13 
20100907 
2013 0915 
2013 10 07 
20140614 
2015 01 18 
201603 14 
2019 0210 
201910 23 
202006 03 
2021 01 03 



56432, 
55186, 
56163, 
56698, 
60137, 
58250, 
58377, 
66131, 
64688, 
60451, 
61639, 
58499, 



674610 
874360 
014200 
641850 
190450 
292000 
929600 
474500 
042880 
414110 
370800 
059467 



2013 05 20 
2009 12 21 
2012 08 23 
201402 10 
2023 07 1 1 
2018 05 11 
2018 09 16 
2039 12 08 
2035 12 26 
202405 20 
2027 08 21 
201901 15 



For WASP-13 and WASP-21 my analysis has resulted in 
the upward revision of the planetary radii, making them two of 
the least dense planets known. For WASP-13 Rh changes from 



1.389 



+0.045 ■ 



, and for WASP- 



XO-2 



21 the change is from 1.143i;j:;;|j; R.jup ^ 1-263 ± 0.085 ± 
0.029 R,Tup- Finally, the inclusion of an improved Tcfr measure- 
ment makes WASP-2 one of the best-understood TEP systems. 

Previously found correlations between Po^b and planet grav- 
ity, and Porb and planet mass, have been revisited and found to be 
of improved statistical significance (3.9(j and 3.8(t respectively). 
The division of TEPs into two classes based on their Safronov num- 
ber is shown to be spurious. 

I recommend further observations of specific types for the ma- 
jority of the objects studied here. These are summarised in Sect.|7] 
In particular, some orbital ephemerides are of low precision and 
will soon be of limited use. The affected systems would benefit 
from new transit photometry to solve this problem. Finally, the dis- 
covery rate and distribution of sky positions of the known TEPs are 
plotted and discussed. 

The main results from this work will be made available in a 
convenient format in the Transiting Extrasolar Planets Catalogue 
(TEPCaQ- 



planet is Pb = 1-73 ± 0.06 Rjup, which is much larger than 
all previous determinations (1.23-1.38 Rjup). This radius makes 
OGLE-TR-56 b one of the largest known planets, as expected for 
its high equilibrium temperature and hot and massive host star. Its 
measured density is smaller by a factor of three. Previous analy- 
ses were based either on survey-quality photometry, a single transit 
light curve which is now known to be strongly affected by cor- 
related noise, or on unnecessary assumptions. Two morals can be 
deduced from this. Firstly, one should always use a newly-derived 
stellar density in calculating revised physical properties of a TEP 
system. Secondly, results based on only one transit light curve can- 
not be definitive as systematic errors may lurk undetected in the 
data. 

New results have also been calculated for OGLE-TR-111, 
OGLE-TR-113, OGLE-TR-132 and 0GLE-TR-L9. These com- 
prise the first complete analyses of recently published high-quality 
photometry for each system. 
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